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We discuss the properties of surface plasmons-polaritons in graphene and describe three 
possible ways of couphng electromagnetic radiation in the terahertz (THz) spectral 
range to this type of surface waves, (i) the attenuated total reflection (ATR) method 
using a prism in the Otto configuration, (ii) graphene micro-ribbon arrays or monolayers 
with modulated conductivity, (iii) a metal stripe on top of the graphene layer, and (iv) 
graphene-based gratings. The text provides a number of original results along with their 
detailed derivation and discussion. 
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1. Introduction 

It is since the days of Arnold Sommerfeld, back to 1899, that plasmonic effects in 
materials and gratings are investigated both theoretically and experimentally. Free 
electron plasma related phenomena in metallic gratings were investigated as early 
as 1902 by Wood ^ . An attempt to a theoretical interpretation of Wood's results 
was first made by Lord Rayleighl^l^ followed, three decades latter, by FanoHH Sys- 
tematic studies of plasmonic effects became possible with the work of Kretschmann 
^ and Otto who devised two different methods, using prisms ^ on top of or be- 
neath a thin metallic film, of exciting surface plasmon-polaritons (SPPs). SPPs are 
evanescent electromagnetic waves coupled to the free electron plasma oscillations, 
propagating along the surface of a conductor. Meanwhile, reliable theoretical meth- 
ods for studying scattering of electromagnetic radiation (ER) by metallic gratings 
have been developed by several authors UMini 

Modern plasmonics gained a renewed interest with the discovery of the anoma- 
lously high transmittance through a periodic array of holes with a size smaller than 
the diffraction limit The effect was explained invoking surface plasmons (SPs) 
|12|13|14|15|16|17[ Another impressive SP-related effect is the surface enhanced Raman 
scattering (SERS) Elllll Beyond the fundamental physics, the plasmonics encom- 
pass a wide range of applications such as spectroscopy and sensing |i8,2i|22|23]24|^ 
photovoltaics 1211, optical tweezers HIZI, nano-photonics IMlll, radiation guiding 1221, 
transformation and Fourier optics IMHEll, etc. These applications rely on the short 
wavelength of SPPs (compared propagating photons of the same energy), strong 
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dependence of their dispersion relation on the environment dielectric constant, and 
high local field intensity associated with localised SPs. 

Doped graphene sustains surface plasmons | 33 | 34 | 35 | 3 6]37] -^j^ose frequency is pro- 
portional to the 1/4 power of the electronic density, a result specific to single-layer 
graphene, and to the 1/2 power of the wave number, a behavior shared with the 2D 
electron gasl^. Moreover, the carrier density in graphene can be varied continuously 
from nearly zero to w lO^'^cm"^ for either type of majority carriers (electrons or 
holes) by the application of an external voltage. This important ambipolar doping 
effect has opened the tantalizing prospect of real time control of SPPs by using 
a gate. In addition, the transfer of graphene films to a range of substrates is rou- 
tine nowadays and has inspired the design of graphene-based structures with unique 
plasmonic signatures, such as quantum dots (or anti-dots) |39 | 40 | 4i | 4 2]43| g^^^^ nanorib- 
■bonsEE313133H In this respect, the question of whether the classical description of 
SPs still holds for finite-size systems, such as quantum dots, is particularly relevant 

m 

The robustness of graphene SPPs with respect to external perturbations, as well 
as the interaction of SPPs with individual quantum systems, are active topics of 
research: the effect of applied stress in the graphene plasmon dispersion E1I121, and 
the role of SPPs in graphene on the decay rate of nano-emitters |5i|52|53|54j55] g^j.g 
now well established. Plasmons in bilayer graphene have also attracted attention 
^ and the transverse electric (TE) mode spectrum has recently been obtained 
The plasmon dispersion relation of a graphene double-layer, two closely separated 
graphene sheets, has also been addressed HHIE A natural extension of the study of 
propagation of SPPs in monolayer graphene is the investigation of the same effect 
in graphene double-layers. In this case, the two degenerate dispersion relations of 
each of the layers hybridize, giving rise to optical and acoustic branches |60|6i]6 2[ 

A series of recent experimental works | 63 | 64 | 65 | 6 6] triggered a revival of interest 
to the plasmonic effects in graphene [MIH In particular, it has been shown that 
graphene has a strong plasmonic response in the THz frequency range at room 
temperature THz photonics is emerging as an active field of research ^ and 
graphene may play a key role in THz metamaterials in the near future. Ju et al. 
have shown ^ that ER impinging on a grid of graphene microribbons can excite 
SPPs in graphene, leading to prominent absorption peaks, whose position can be 
tuned by doping. Fei et al. mapped the plasmon dispersion relation of graphene, 
and demonstrated the plasmonic tunability of graphene In agreement with the 
theory,!^ infrared SPPs in graphene have been shown to possess remarkably large 
propagation lengths when gauged against more conventional structures A simi- 
lar experiment was performed by Chen et al. where excitation and subsequent 
detection of SPPs were achieved. Yan et aL'^have shown that graphene/insulator 
stacks can be used as tunable infrared plasmonic devices, able to work both as a 
filter and as a polarizer, having the potential for far-infrared (FIR) and THz pho- 
tonic devices. As a proof of principle experiment, it has been shown that this type of 
devices can shield 97.5% of ER at frequencies below 1.2 THz^^. Last, the possibility 
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of transforming graphene into a mantle cloak, working in the THz spectral range, 
has also been addressed'^. 

In general, SPPs cannot be excited by directly shining light on a homogeneous 
system due to kinematic reasons: the momentum of a surface polariton is much 
larger than that of the incoming light having the same frequency. Therefore, some 
type of mechanism is necessary to promote the excitation of SPPs. The most com- 
mon mechanisms for SPPs excitation are: (i) attenuated total reflection (ATR) 
(ii) scattering from a topological defect at the conductor surface IHIZIl, and (iii) 
Bragg scattering using diffraction gratings 1^ or a periodic corrugation of the sur- 
face of the conductor HEH The method of Ju et al. is similar (but not identical) to 
the patterning of a metallic grating on top of graphene. A theoretical account of 
the experiment by Ju et al. was given in a recent workS^l. A novel method of SPP 
excitation consists in using the apex of an illuminated nanoscale tip this method 
provides a two-orders of magnitude enhancement of the in-plane momentum rela- 
tively to that of the impinging ER in free space bypassing the referred kinematic 
limitation. It has also been shown that modulation of the optical conductivity gives 
rise to efficient ER coupling to SPPs in graphene without the need of a grating. 
The same principle works with split ^ or modulated 1^ gates. Free space excitation 
of SPPs using non-linear optics effects has been proposed EEll and effectively put 
in action 1221. 

The reason why a periodic corrugation allows for the excitation of SPPs can be 
understood in analogy with the theory of electrons in a periodic potential. Here the 
periodic corrugation plays for SPPs the same role as the periodic atomic potential 
plays for electrons, that is, the SPP momentum is conserved up to a reciprocal 
lattice vector. The periodic corrugation provides the missing momentum needed to 
excite the polariton. Another way of understanding the effect is to note that the 
grating gives rise to a SPP band structure. Then the folding of the SPP dispersion 
curve, forming bands in the first Brillouin zone, makes it possible for an incident 
electromagnetic wave to excite a SPP mode associated with the upper bands for the 
same wavevector. The excitation of SPP modes in the first band still is not possible 
in the grating configuration unless the ATR technique is used. 

This paper is organized as follows. In Sec. [2] we write down the macroscopic 
form of Maxwell's equations, introduce boundary conditions at an interface con- 
taining a graphene sheet, and provide the definition of TM (transverse magnetic) 
and TE waves. In Sec. [3] we review the basic properties of the optical conductivity 
of graphene. In Sec. |4] we compute the spectrum of a TM SPP wave providing both 
numerical and analytical results. In Sec. [5] we derive the dispersion relation of a 
TE SPP wave. In Sec. [6] we discuss in detail the excitation of SPPs on graphene 
using a prism in the ATR configuration originally proposed by Otto In Sec. [s] 
we discuss scattering of ER by an array of graphene microribbons and other fiat 
graphene structures with periodically modulated conductivity. Section [7] is devoted 
to the SPP excitation by shining light on a metal stripe deposited on graphene. In 
Sec. [9] we give a formulation of Rayleigh approximation to scattering of ER by a 
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graphene-based grating due to Toigo et al. adapted to this particular geometry. 
There the problems of both a sine a sawtooth gratings are solved. In Sec. [lO] we 
discuss SPPs on a metal grating coated with graphene. In Sec. [iT] we summarize 
the main results of the paper. In all sections we give enough details, so that the 
interest reader may reproduce all the results by him/herself. 



2. Maxwell's equations and boundary conditions 
2.1. Planar geometry 

The central idea of the emerging area of nanoplasmonics is the use of SPPs, coupled 
charge-radiation excitations existing at a dielectric/conductor interface, for appli- 
cations ranging from chemical sensors and surface-enhanced Raman spectroscopy 
(able to detect a single molecule) to solar cell's optimization. One of the most at- 
tractive features of SPPs is that they concentrate and guide ER at subwavelength 
scales. This is appealing since one can conceive circuitry using metals embedded 
in dielectrics, which is able to propagate both electric signals and SPPs In 
graphene, the SPP wavelength can be about 40 times shorter than the wavelength 
of the impinging radiation in free space 
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Fig. 1. Tangential and perpendicular fields at an interface. A graphene sheet is located between 
the two dielectrics. 



We shall consider a graphene sheet, with conductivity a, cladded by two di- 
electrics of relative dielectric permittivity ei and £2- Graphene can be seen as the 
ultimate conductive thin film able to support SPPs As in any other conductor, 
the free charges in graphene couple to ER. We seek solutions of Maxwell's equations 
in the form of surface waves propagating along the graphene sheet. As we will see, 
graphene supports two different types of surface waves, the TM (or p— polarized) 
and TE (or s— polarized) SPPs. 
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2.2. Maxwell's equations 

In the MKS system of units, the macroscopic Maxwell's equations read 

V-D^Pf, (1) 
- - dB 

Vxi.^--, (2) 

V • B = , (3) 

^ ^ ^ dD 

VxB^J^ + — , (4) 

where D = eegE, B — iiqH , e is the relative dielectric permittivity, cq and /io are 
the vacuum dielectric and magnetic permeability, p/ is the charge density (charge 
per unit volume), and J/ is the current density (current per unit area). 

In the case of a 2D metal, such as graphene, and assuming that it lies in the 
xy— plane, we have 

Jf = JsS{z) , pf = psS{z) , (5) 

where Jg is the surface current density (current per unit length), ps is the surface 
charge density (charge per unit area). For a 2D metal with a frequency-dependent 
(optical) conductivity cr, we have Jg = aEt, where Et lies in the xy— plane as well. 

The boundary conditions at the interface between the two media are (see also 
Fig.[l]): 

Eit = E2t , (6) 
Hit -H2t^-fs-xn, (7) 
Dm - D2n = Ps , (8) 
Bm - B2n ■ (9) 

Since graphene is a two-dimensional system, it enters in the calculation of the surface 
wave dispersion relation only through the boundary conditions. Then, the only 
quantity we need to known is its optical conductivity. 

Assuming the time dependence of the fields in the form e~*"* and considering 
the absence of free volume currents and charges, Eqs. ^ and Q are written as 

dyE, ~~ d,Ey - ILUB, , (10) 

dzEx - dxE:, = iujBy , (11) 

d.,Ey - dyE, = iojB, , (12) 

and 

dyBz — dzBy = —ic~'^eujEx , (13) 

dA - d,B, = -ic-\LuEy , (14) 

d^By - dyB^ = -ic^'^euEz , (15) 
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respectively. As mentioned above, Maxwell's equations apply to the dielectrics with 
permittivity ei and €2 surrounding graphene, and we take into account the relation 



2.3. Definition of TE and TM modes and Poynting vector 

A conductive surface, depending on certain properties of metal's optical conductiv- 
ity, can support transverse electric (TE or s— polarized) or/and transverse magnetic 
(TM or polarized) surface waves. These waves propagate along the metallic in- 
terface and decay exponentially away from it (along the z— direction). In the case 
of graphene, the free carrier oscillations are confined to the xy plane, while the 
electromagnetic field penetrates considerably into the cladding materials. 
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Fig. 2. Electric and magnetic fields of transverse electric (TE; left in the figure) and transverse 
magnetic (TM; right in the figure) surface waves. The graphene sheet lies in the xy— plane. 



In Fig. [2] we represent both types of waves propagating along the x direction 
with momentum q. For TE- waves the electric field is oriented along the y direction. 
Then, we have 

X (TyyEyij X fl^ <7yyEyX , (16) 

and for TM-waves the magnetic field is oriented along the y direction, so we have 

X (TxxExX X h = -cFxxExy , (17) 

where we have assumed an anisotropic response. For an isotropic system, such as 
unstrained graphene, axx = o'yy = a. 

For future use, we recall here the definition of the Poynting vector, which is the 
instantaneous ER energy flux per unit area 

S^ExH, (18) 

where E and H are assumed to be real. Representing E and H by harmonic functions 
we have: 

S ^ ^{Eoe''''-"^*) X 5R(i?oe"'^~'"*) 

= ^^{Eo X H*) + ^^{Eo X H^e^^'^^^-^'^) . (19) 
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The time average of S, i.e. the directional energy flux density, reads: 

{S)=^-^iEoxH*). (20) 

The calculation of (5) allows for the determination of the ER reflectance and trans- 
mittance at an interface. 



3. Optical conductivity of graphene 
3.1. General expression 

In this section we provide a brief overview of the optical properties of graphene. 
Its optical conductivity is a sum of two contributions: (i) a term describing in- 
terband transitions and (ii) a Drude contribution, describing intraband processes. 
At zero temperature the optical conductivity has a simple analytical expression 
EMZMEMSIH The inter-band contribution has the form cr/ = cr^ + ia'^ , with 

, / 1 hLj-2EF 1 huj + 2EF\ , , 

= do 1 H — arctan arctan , (21) 

and 

1 (2g/. + M^ + ^V 

a, = -do - In ^^EF-hur + h^j^ ' ^^^^ 

where do — 7re^/(2/i), e < is the electron charge, 7 is the relaxation rate, and 
Ep > denotes the (local) Fermi level position with respect to the Dirac point. 
The Drude conductivity term is 

o-n=do — . (23) 



The total conductivity is therefore given by 



ao . (24) 



In Fig. |3] the two contributions, Drude and interband, are plotted separately for a 
given value of Ep and F = ^17. For heavily doped graphene and for photon energies 
Huj/Ep <C 1, the optical response is dominated by the Drude term. Therefore, in 
what follows we assume 

dg « dD , (25) 

since we are interested in this regime of frequencies. For the frequency range of 
interest in this work (the THz spectral range) the above approximation gives accu- 
rate results. The only exception to this approximation will be made when discussing 
TE- waves in graphene (Sec. Isl). 
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Ep=0.45 eV; r=2.6 meV Ep=0.45 eV; r=2.6 meV 
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Fig. 3. Optical conductivity of uniform graphene: Drude (left) and inter-band (right) contribu- 
tions. We assume Ep = 0.45 eV and F = 2.6 meV. The solid (dashed) line stands for the real 
(imaginary) part of the conductivity. In both panels uip = Ep/h. 

3.2. Drude conductivity in a magnetic field 

When a static external magnetic field is considered, the response of two-dimensional 
electronic systems to an external time-dependent field is described by the magneto- 
optical conductivity tensor aap, where a,l3 = x,y denote the in-plane coordinates. 
The tensorial nature of the magneto-optical response is a direct manifestation of 
the Lorentz force which, in the presence of a magnetic field, gives origin to both 
longitudinal and transverse electronic currents. 

We assume that the applied static magnetic field is homogeneous and transverse 
to the graphene sheet, B = Be^- In the semiclassical regime, huj/Ep <^ 1, the 
magneto-optical transport in graphene can be described in terms of Boltzmann's 
transport equation. ^ Within this formalism, the electric current is obtained ac- 
cording to 

J'=^ J d'kSfj:Vf:, (26) 

where is the deviation of the carriers' (electrons or holes) distribution function 
from the equilibrium Fermi-Dirac distribution, /o(e), and 

Vj:^ VF{cos9j:,smej:) , (27) 



with vp « 10^ m/s denoting the Fermi velocity in graphene and 0g = arctan(A:j,/fc3;). 
We remark that both spin and valley degeneracies have been included in Eq. (26). 
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In the presence of an electromagnetic field, the distribution function perturba- 
tion, (5/g, is the solution of the kinetic equation ISH^ 



where the standard relaxation approximation, 

(29) 



dt 



scatt 



has been assumed. Here denotes the carrier's relaxation time, and E^^ — {Ex,Ey) 
is the projection of the electric field onto the graphene plane. The kinetic equation 
can be solved exactly by writing as 

<5/,, = e— (30) 

and noting that (wg x B) ■ dj^dfj: — ■ {B x d^6J-j^), in order to obtain after some 
algebra 

1 /I - iujTt -TtUJr. 



^'^ (1 - iiUT^f + LoItI \ T^,OJ, ' 1 - UtJ ' ^^^^ 



where £ is defined as 



and 



£^i:-~^^iE.vj:^.'EyV^,.,y), (32) 



uJc = evlB/\EF\ (33) 



is the cyclotron frequency in graphene. Introducing the explicit form of Sf^ in 

S{e ~ Ef) , (34) 



Eq. (261, and assuming T = 0, 

_dfo 
de 



Cxx{B,uj) = . , , ,^,^2 , (35) 



a^y{B,uj) = - — — ———2, (36) 



one arrives at the semiclassical form of the conductivity tensor: 

e^2\EF\ "i~iuj 
'h h {^-iuY + [Lo^[B)y 
2Ef 0Jc{B) 
h h (^^ iuj)2 + [uj^{B)] 

Uyy(B,ui) = Uxx(B,uj), and ayx{B,uj) = —axy{B,ui). Note that for simplicity we 
have expressed the result in terms of the relaxation's rate 7 = . 

In the presence of a magnetic field, and for pristine graphene (/17 energy 
scales), the Drude peak (i.e., the maximum of YUia^x) is located at the cyclotronic 
frequency, reflecting the intuitive fact that the optical response (absorption) is max- 
imum for impinging light in resonance with the frequency for cyclotronic motion, 
that is, uj = LJc- We note that the Drude conductivity for zero-field [Eq. (23l] is 
recovered by setting i3 = in Eq. (35). 
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4. Spectrum of TM SPPs in graphene 
4.1. Dispersion relation 

Let us find the form of a p-polarized surface wave in graphene. We assume a solution 
of Maxweh's equations in the form 

E={E„,^,,0,E„,,,)e'''-e-^^\''\ , (37) 
5 = (0,S,„,^,0)e^«^e''="l^l, (38) 

where m = 1, 2 refers to the media 1 and 2 (see Fig. [ij, waves are exponentially 
decaying in both directions away from the graphene sheet, and we seek the dispersion 
relation, uj = uj{q), of this type of waves, with q denoting the wave number along 
the graphene sheet (see Fig. [2]). The wave numbers yet to be determined. In 



this case, Maxwell's equations (11), (13 1, (15) yield 



KmKi^x - iqEm^z = i^B„i^y , (39) 

(-l)"+V„,B„,j, = -iuJC-^€^E„^^^ , (40) 

qB„i^y — ~ujc~^emE„i^z , (41) 

which can be solved in terms of the amplitude of the magnetic field and allows for 
the determination of k,„, that is, 

i^™.. = * — i?™,,(-l)™+\ (42) 

UJem 

Em.z = —-^—Bm,y , (43) 

4, = q^- u\„Jc^ , (44) 

with B,n y a constant. The SPP spectrum follows from the boundary conditions (|6]) 
and 0: 

£'1,2; — E2^x <^ Bi y — —B2^y , (45) 

Bi^y — i?2„ — o-xxEi^x , (46) 

from which we obtain the dispersion relation (in an implicit form, since both 
and depend on frequency) 

1 + — + idxx z — 77- = ; 

. (47) 
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We note that Eq. (47) has real solutions only when the imaginary part of the 
conductivity is positive. This takes place when the conductivity is dominated by 
the Drude contribution. If the real part of the conductivity is finite (non-zero), the 



solutions are necessarily complex. Eq. (47) gives the spectrum of the polarized 
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SPPs in graphene. When the two media are the same (ei = £2 = e), we obtain a 
simpler relation for the spectrum, 

l+i_^^^^g2_^2g/c2 =0. (48) 

Notice that starting from the boundary condition (|8| would lead to the same 
spectrum. It can be shown as follows. From the continuity equation 

dtPsix,t) + V ■J,{x,t)=0, (49) 

follows —iujpg + iqUxxEi^x = 0- Thus, the boundary condition reduces to Ei z^-i — 



E2,z^2 — f^x'^xxEi^x/i^^o)- Using the relation (431, we obtain 



— Bl^y + i?2,y — (7xxEix , (50) 



which is the same as Eq. ( 46 1 



4.2. Simplified analytical form 

Let us obtain some simple analytical results for the spectrum of the SPPs in 
graphene. We assume that the conductivity of graphene is given by the Drude 
contribution only. Ignoring absorption (7 = 0), we have 

ao^i-: V = (Jo — ;- • (51) 



Furthermore, assuming £2 = ei = e, the equation for the SPPs spectrum reads 

1 1/ 



(52) 



Inverting and squaring, we obtain 

„4 ^ f J^Vr„2_,„2^, 2^ 



When cj — > 0, we obtain from Eq. (53 1 uj — vq, where v = c/y/e is the speed of light 
in the dielectric. In the electrostatic limit (also dubbed non-retarded or plasmon 
approximation) we consider uP' jv^ <^ q^, in which case we obtain 

2 2 2aEF . 

h^u}'^ = q- = hcq , (54) 

2eeo e 

where 

1 

(55) 



Aneohc 137 



is the fine-structure constant. We note that Eq. (54 1 coincides with the spectrum of 
plasmons in graphene. Since Ep = vphkp, the spectrum depends on the electronic 
density as 71^^, which is specific of graphene and was experimentally confirmed 
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We now derive the spectrum of the SPPs taking into account absorption. In this 
case we have to use the fuh form of cr^i. The dispersion relation now reads, 

K = -^(w^ + i7cj) . (56) 

In the electrostatic limit we write k~ q = q' + iq" (assuming q" <C q')^ which is now 
a complex quantity due to a non-zero 7. In this case, the SPP spectrum is obtained 
from 

q' + iq" ='^-^{uj^ + l^Uj), (57) 

V 

where we have assumed 7^0;. In terms of q', the spectrum is 

(58) 



which has the same functional form as in Eq. (54). The decay of the wave as it 
propagates in space is characterized by g", 

q" = ^7^, (59) 

from which follows 

(60) 

q Lo 

that is, the decay is less pronounced at higher frequencies and smaller 7. The con- 
finement of the electromagnetic field in the z direction (i.e. its penetration depth in 
the dielectrics surrounding the graphene sheet) is given by . In the electrostatic 



limit, from Eq. (54) we have: 

K ^ q = u)^ ^ . (61) 



2 2eeo 



In the dielectric the wave number of light is k — lj/v. Comparing k, with k we 
obtain: 

K 2eeQV Je hiu 137 ^luu 

— uj — ve • 62 

Thus, the EM penetration depth decreases with e and with the photon frequency 
and increases with the Fermi energy. Considering, as an example, w/(27r) = 10 
THz, Ef — 0.2 eV, and e = 4, we obtain k^^ « 0.035fc^^, a fairly high degree of 
confinement in comparison with the wavelemgth. 

4.3. Numerical results 



In the general case of different dielectrics, neglecting absorption, Eq. (47) can be 
written as 

ei 62 AoEf 



(63) 
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q(1/nrn) q ( 1/|^m ) 



Fig. 4. SPP dispersion curves calculated taking into account only the imaginary part of the 
conductivity (thus q is real). We have chosen h'y = V = and Ep = 0.45 eV. Left: ei = 3 and 
£2 = 4; right; ei = 1 and t2 = 11 (silicon). The straight (dashed) line stands for the light dispersion 
fiio = hcq/,Je2 in the medium 2 (that with higher dielectric constant). 

which is convenient for numerical purposes. [In the numerical calculations we take 
he = 0.2 eV-/zm; it is also useful to recall that w/(27r) =1 THz corresponds to an 
energy of 4.1 meV and to a wavelength of 300 /J,m.] 



In the non-retarded approximation it is possible to derive from Eq. ( 63 1 analyt- 
ical results which coincide with those of the previous section upon the replacement 
e — >■ (ei -I- £2)72 = e, that is, 

hVLp K y/2aEFhcq/l. (64) 

Note that we have used the symbol 57p = ilp(q) to denote the explicit solution of the 
dispersion relations; we shall keep this notation throughout the paper. In general, 
Eq. (|63[ ) has no analytical solution. In Fig. |4] we give the numerical solution to 
Eq. ( |63| in two different cases. 

In Figs. [5] and |6] we present the low wavenumber part of the SPP spectrum. 



comparing the numerical solutions of Eq. (47 1 obtained with different levels of ap- 
proximation to the optical conductivity and considering the role of the damping 
(or homogeneous broadening) parameter, F = fvy. In Fig. [s] the effect of the damp- 
ing is demonstrated, taking into account the full conductivity ud- In this case the 
wavenumber q is complex, q = q' + iq" , where g" describes the decay of the SPP 
as it propagates in space along the graphene sheet. From this figure we see that 
the effect of the increase of the damping is two-fold: it shifts the SPP dispersion 
relation toward higher energies, for the same q' , and enhances the value of q" (as 
expected). If one wants to have long propagation lengths for the SPPs, then F must 
be as small as possible. In the same figure we also represent the light-lines hcq/y/e[ 
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Ep=0.45eV, e1=2, e2=3 




q' (1/nm) 



Fig. 5. Plasmon-polariton dispersion curves for ei = 3, e2 = 2, considering the effect of broaden- 
ing (for different values of F). Full expression for the conductivity with Ep = 0.45 eV was included 
in the calculation. The red dashed line is for zero broadening. The inset shows the dependence of 
the imaginary part of q (g") on the real part of q (q') (notice the logarithmic scale of the vertical 
axis). 



and hcq/y/e2, which correspond to the dispersion relations for photons propagating 
in the media 1 and 2, respectively. It may seem that one should be able to excite 
SPPs in graphene by directly shining ER on it if F is sufficiently large because 
the dispersion curves in Fig. [5] intersect the light lines, as most clearly seen for the 
curve corresponding to F = 26 meV. This is not true. For such a high value of F, 
plasmon-polaritons are overdamped (note that q' < q"), the "dispersion curve" uj 
versus q' — ^{q) is not quite meaningful, and the intersections do not have their 
usual physical meaning. 

In Fig. [6] we present the effect of neglecting the real part of the optical conduc- 
tivity. 

The central feature is the vanishing of the dispersion curve for a finite value of g', 
which can be expressed as 

, £1+62 h'^^ 

Q = -r-F^ ■ (66) 

This is a spurious result. When we use the full Drude conductivity expression, the 
dispersion relation vanishes at zero q' . Changing the dielectric constants of the 
surrounding media changes the value of q' for the same frequency, as can be seen 
comparing the left and right panels of Fig. [6j 



4.4. Structure with two graphene layers 

In the previous sections we have studied the SPP spectrum of single layer graphene. 
We have shown that the typical value of SPP's energy for wavevector q ^0.1 /im~^ 
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r=2.6 meV; Ep-0.45 eV; ei=i, e2=l r=2.6 meV; Ep-0.45 eV; el=2. e2=3 




q'(1/Mm) q'(1/Mm) 

Fig. 6. SPP dispersion curves calculated considering the real part of the conductivity either 
finite or zero (but keeping 7 finite in the imaginary part) for cases of equal (left: ei = 1, e2 = 1) 
and different (right: ei = 3, £2 = 2) dielectrics surrounding graphene with chemical potential 
Ep = 0.45 eV and damping h'y = 2.6 meV. The straight lines marked ei and £2 refer to the light 
dispersion huj = hcq/^/ei and hu} = hcq/^/e2 in the dielectrics 1 and 2, respectively. 

is huj ~4 meV, corresponding to a frequency of about 1 THz. We would like to 
investigate the possibility of shifting the resonance frequency towards higher values. 
One way of achieving this is using a double-layer structure as shown in Fig. [Tj The 
two degenerate SPP branches, associated with each of the layers, hybridize giving 
rise to a two branch spectrum, with one branch having higher energy than the bare 
spectrum for each layer. 




Fig. 7. Structure with two graphene layers separated by a dielectric with £2- Field profile for the 
"acoustic" SPP mode is shown qualitatively. 

We first compute the spectrum of the double layer system. We assume a configu- 
ration as represented in Fig.[7]with semi-infinite media 1 and 3 and a dielectric layer 
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2 of thickness D. Applying the boundary conditions at two interfaces and requiring 
exponential decay in the infinity yields: 



K3 



-1 


1 


1 


+ i — 


_ 12_ 


12. 


ujeo 


K2 


K2 









K2 


ft2 





gK2-D 





+ i 







E3,x 












E2-,x 









(67) 



The spectrum is given by the vanishing of the determinant of the above matrix, 
that is, 





(7 

f i 






l-i 








K2/ 






'«2 



-K2D I 



+ i h 2 

Ki weo K2 / V '*3 t^eo '«2 




Fig. 8. Spectrum of the structure of Fig. [t] with its lower (acoustic) and upper (optical) branches. 
The dashed curve is the spectrum of a single layer. The parameters are em = 1 (with m = 1, 2, 3), 
D = 1 fira, and Ep = 0.45 eV. 



As expected, in the limit q2d — > 00, the two interfaces decouple and Eq. (68 1 
reads as 





a 

f I 


^'1 


r- 


cr 

f i 






ujeo 


H2J 


\l^3 




K2 



= 0. 



(69) 



Each term in brackets represents the plasmons-polaritons spectrum of one graphene 
layer [c.f. Eq. (47)], degenerate if ei = £3. 
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In the "symmetric" case of ei = £3, the dispersion relation (68) can be factorized 



cos\i{k2D/2) 



£2^ 
K2 



sinh(K2-D/2) 



\{-- 









K2 



(70) 



The first product term in Eq. ( 70 1 corresponds to the symmetric mode with the 
electric field z dependence even with respect to z = D/2, while the second term 
describes the antisymmetric mode, odd with respect to the central plane. 

Another analytic limit is that of ei = £2 = £3 = e and q ^ uje^^^ /c (non-retarded 
or plasnion approximation) 



In this case, Eq. (|70p reduces to 
2ee'^ 



qa 



0, 



(71) 



where upper and lower signs correspond to the first (symmetric) and second (an- 



tisymmetric) terms in ( 70 ) , respectively. If we assume that the conductivity of 



graphene is given by Eq. (51 1, it is possible to obtain a close form for the spec- 
trum of the surface waves. 



2o2 



2aEF 



hqc (1 ± e 



-qD\ 



(72) 



Note that for qD > the spectrum given by Eq. (72) becomes inaccurate since 
the condition qD ^ wge^/^/c is violated. In this regime the spectrum can only be 
obtained numerically. In Fig. [s] we give the numerical solution of Eq. ( 68 ) and obtain 
the two branches of the spectrum; clearly they both vanish at zero q. The optical 
(symmetric) branch lies above the dispersion curve of a single graphene layer and 
the acoustic branch is below it. This dependence can be obtained analytically as 



well by inspection of Eqs. (72) and (54) 



The propagation of TM-waves in multilayer graphene-based structures, under 
strong light illumination, was considered by Dubinov et al. 



5. Spectrum of TE SPPs in graphene 

A broadband polarizer based of on the propagation of TE-surface waves on graphene 
has been demonstrated It is, therefore, interesting to compute the properties of 
a TE surface-wave. Let us find the spectrum of an s— polarized (TE) surface wave 
in graphene. These type of waves do not exist in the traditional 2D electron gas 
because the imaginary part of the conductivity is always positive. We assume a 
solution of Maxwell's equations in the form 

B = (B„,,„0,i3„,,)e'«^e-''™l^l, (73) 
^= (0,i?„,„0)e^«-e-''"l^l, (74) 
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where m = 1, 2 refers to the media 1 and 2 (see Fig. [ij. In the case under consid- 
eration the boundary conditions ^ and ([?]) are represented as 

El^y — E2,y , (75) 
Bi^x — B2,x = f^Of^yyEi^y ■ (76) 

For this particular case, Maxwell's equations read 

( 1) ^ l^mEm.y = i^B^ x ; 

iqEm,y = iujBm,z , 
(~1) l^mBm,x ~ i^B^i^ ;^ = —ic Cjn^Ejy^ y , (77) 

and we obrain the relations: 

B„,,x^i~ir+'^ — Em,y, (78) 
UJ 

Bm,z — EjYi y , (79) 
up- 

'^m =9^ T^rn ; TO = 1, 2 . (80) 

From the first boundary condition it follows E\^y = Ei^y and the second boundary 
condition ( 76 1 yields the spectrum of the TE- waves 



Ki + K2 - ioj^oayy = . (81) 



We note that Eq. (81) has real solutions if and only if the imaginary part of the 
conductivity is negative (and the real part of the conductivity is zero) . For cr' > the 
wave is damped by the imaginary part of either q or lo, depending on the excitation 
conditions. Since the imaginary part of the Drude conductivity is positive, this 
contribution alone cannot give rise to TE-waves. On the other hand, the interband 
contribution to ct" is negative. Therefore, when the two contributions, for a given 
frequency, add to a negative number, TE SPPs can propagate in graphene. 



Eq. (81 1 cannot be solved analytically. However, the third term is small because 
it is proportional to the fine structure constant. Then, in the simple case where 
ei = 62 = e, the spectrum of the TE wave is essentially equal to huj < hcq/y'^., that 
is, the dispersion relation of a free wave in a dielectric. Only close to the threshold 
for interband transitions, hoj = 2Ep, does the spectrum deviate considerably from 
this result. 



6. Excitation of SPPs by evanescent waves 
6.1. ATR configuration 

As discussed above, excitation of SPPs on a fiat metal surface is not possible by 
direct illumination (excepting the case of a periodically modulated conductivity; see 



Sec. 8.2). It is because the wave vector of the SPP, at a given frequency, is much 



larger than that of the impinging radiation. One way out is excitation by evanescent 
waves. We assume an attenuated total reflection (ATR) conflguration of the form 
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Fig. 9. Schematic representation of the experimental ATR setup needed to excite surface plasmon- 
polaritons in graphene (Otto configuration). The graphene layer is located between two dielectric 
media of relative permittivity ei (considered semi-infinite) and ei (of thickness d). On top of the 
latter there is a prism of relative permittivity £3 > (ei, £2) (which is a necessary condition for total 
internal reflection). The incident angle of the incoming p— polarized wave is Q, and the electric and 
magnetic fields are E = (E^, 0, E^) and B = (0, By, 0), respectively. 



depicted in Fig. [9] Electromagnetic radiation coming from the prism impinges into 
the interface 3-2 at an angle 9 larger than the critical angle for total internal reflec- 
tion, 9 > arcsin Y^max(ei, e2)/e3 ■ So, only evanescent waves can exist in the layer 
with £2 (of thickness d) and in the half-space ei. In this configuration it becomes 
possible to couple the incident wave to the SPPs in graphene because the dielectric 
constant of the prism is larger than those of the dielectrics cladding graphene. 



6.2. Fields in the ATR regime 



We start by deriving the fields in the media 1, 2 and 3 (refer to Fig. 10). The 
incoming wave vector is fc, = {k sin 9, 0, k cos 9) and the reflected one at the interface 
is kr — {ksm9,0,~kcos9). Due to translational invariance, the component q = 
k sin 9 is conserved at all interfaces. 

We study the case of TM-waves, that is, the fields have the form B„i = 
(0, -B„i,i/, 0) and Em — {Em^x,^, E,-a,z) = 1 — 3). If wc further assume that 
i?3 = i3(^'-)e*'^' --^ and £^3 = ij('^'')e'^'.'-*', it follows from the Maxweh's equations 
that 

fc^weg^Vc, (82) 



= ±i^£;(^'-) , (83) 
^ c^k 



^ ^A^C^-) . (84) 
kz 



Here k^ — kcos9. In Eqs. (83)-(84| the upper and lower signs correspond to the 
incident and reflected waves, respectively. On the other hand, if we assume that 
Bm — -B^^e'^^e^'^"'"^ and Ern = E^^ e'-'^^e'^'^'"^ , a situation corresponding to 
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evanescent or exponentially growing waves [simple generalization of Eqs. (37)- (38 1] 
we find 



Si 



±i- 



and the dispersion relation is equivalent to (44). Since q is conserved, we have: 

, ,2 , ,2 \ 1/2 

> Q <^ s.in6 > 



^2 = sm 



£3 



(85) 
(86) 

(87) 



The last condition means that evanescent waves can occur in the medium 2 if and 
only if £3 > €2 - For the geometry of Fig.lTO) the boundary conditions at the interfaces 



graphene 




Fig. 10. Three dielectrics with £3 > ei,e2- Between the dielectrics 1 and 2 there is a graphene 
sheet. 



z = d and z — read as 



E„ 

Br, 



rn+l.x : 



I Em ,X ; 



(89) 



respectively, where 5m,m' is the Kronecker symbol. 



6.3. Total reflection 

Let us assume that d ^ 00 and the graphene sheet and the medium 1 have no effect 
on the reflection at the interface 2-3. We shall show that, for 9c = arcsin (£2/63)^^^, 
the reflection is total, even though there is an evanescent wave in the medium 2. For 
this particular case the boundary conditions are written as Ex + Ex ~ E2 J and 
By^ + By^^ = B^^'y , where the subscripts i and r stand for incident and reflected 
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waves as before. Explicitly, the second boundary condition reads 

- ^4'-) = -^'^Ei-J , (90) 



with the following solution, 

e"'" , (91) 



where we have defined a = arctan(e2fcz/e3K2). Note that for grazing incidence a = 0, 
while for 9 = 9c we have K2 = and a = tt/2. We thus conclude that total reflection 
(|rp = 1) occurs for 9 > 9c- The amplitude of the evanescent field reads 

Ei^^ <^3K2 ~ ie2kz e3H2 - 1^2^ cos 9 ' 



where k,2 — sin 6* — e2/e3 < k. We show in the next section that the presence of 
a graphene layer, as in Fig. 10 can frustrate the total reflection in a dramatic way. 



6.4. Otto configuration: p— polarized wave 
6.4.1. Reflection coefficient 

We shall now compute the reflection coefficient in the situation represented in 
Fig. [TOj where we have three dielectrics, a graphene layer at the 2-3 interface, and 
9 > 6c- We will see that in this case |rp can be smaller than unity (down to zero) 
due to the excitation of SPPs waves in the graphene layer. 

In the dielectric 3 we have both incoming and reflected waves, while in the di- 
electric 2 we can have both evanescent and exponentially growing solutions because 
this region has a finite thickness d. In the dielectric 1 only evanescent waves can 
occur. 

The fields in the three regions are: 
Region 3 

B3,y = B^e^^-^ + B^^e^^-\ (93) 
^3 - 0, e^^^-^+ 0, E''^^^ e''""-^. (94) 

Region 2 

B2.y = B^+^e^^^^'-"" + B^-)e'^^"'•^^ (95) 
E2 = {Ei%\0,Ei^^) e^'-'+ ^'-'^ (96) 

Region 1 

5i,.-si7;e^^T'-, (97) 
i?i = (<j,0,i?(-')e^^l"-'\ (98) 
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where km,± — {ksm6,0,^iKr, 
z = and z — d are: 



1,2. The boundary conditions at the interfaces 



?(+) 



y ^ y 



B 



2,!; 



■2,V ■ 



(99) 
(100) 



(101) 
(102) 



ExpUcitly, using relations (83) and (85), we can rewrite Eqs. (100) and (102) in 



terms of the x-component of the electric field, 
2 = 0: 

h ^ k ^ 

z = d: 



1^2 



1^2 



K2 ■ K2 



a 

+ 1 — 



-Kid 77i(-) 



(103) 



(104) 



After the elimination of e[ J , the linear system of equations defined by the boundary 
conditions can be written as 

'Ei^^/Ei^y 

E2,x I Ex 
.E2,x/Ex 



P' 



(105) 



where 



-1 
K2e3 



1 

ik^e2 



1 

-ikze2 



e«^'^(K2A^+^ + xi) e-«^'^(K2A^+^ - xi) 



K2e3 




A(±) 



(we™ ± iaK^/eo) ■ 
The solution of the linear system yields: 

K2e3??i + ikzf-2ri2 



K2e3f?i - ikz^2ri2 ' 



Kid TT"^ 

£jx -Cjx 



E. 



(-) 



P(0 



2cje2KiK2e3 
K2e3?7i - ikze2m 



(106) 



(107) 

(108) 
(109) 



(110) 
(111) 



where 771 = '«2Ai^^ sinh(K2rf) + cosh(K2d)xi and 772 = K2AJ'''^ cosh(K2(^) + 
sinh(K2(i)xi- We note that, if the conductivity has an imaginary part only, then 
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the function ' is real and the reflection coefhcient is simply a phase. Therefore, 
the reflectance is equal to unity. Thus, for detecting the excitation of a surface wave 
the conductivity of graphene must have a finite real part. In Fig. [TT] we represent 
the reflectance and the absorbance as functions of the energy of the incoming ER. 
A sharp dip in the reflectance spectrum is seen at an energy of about huj 4 meV. 




(0 (meV) 



Fig. 11. Absorbance and reflectance spectra of graphene in the Otto configuration. The param- 
eters are; Ep = 0.45 eV, V = 0.1 meV, ei=14, e2 = I, e-j, = 5, d = 2.51 /^m, 6»c = 36.7°, and 
6i = 2.256c. Sohd black hne: |rp; Red dashed hne: l-|rp. In the calculation we have used the 
Drude formula for the conductivity of graphene, Eq. \23\. 



We also note a reduction of the reflectance close to zero energy. We want to un- 
derstand why the reflectance falls to zero at these two energy values. First, we make 
a three-dimensional representation (Fig. 12), where the absorbance is represented 
as function of both the angle of incidence and the incoming photon energy for a 
given width d. Maximal absorbance (minimal reflectance) occurs close to zero and 
4 meV, and is accompanied by a drastic change (from positive to negative values) 
of the phase of the reflection coefficient as depicted in Fig. [T2j The comparison of 
top and bottom panels in Fig. [12] shows that the resonance maximum of absorbance 
corresponds to maximal amplitude for excited SPPs and phase ss — 7r/2. The latter 
imposes the fulfillment of the condition Ke{E[~J e~'^^'^) w at resonance, which, 
after taking into account Eq. (Ill), leads to the following constraint: 

Re(?7i) = I — - — I sinh(K2d) + — cosh{K2d) = . (112) 

\Kl UJeoJ K.2 



Eq. ( 112 ) determines the frequency onset Wniin' main idea can be interpreted 

in the following manner. In the limit d — > cx) this equation can be obtained by the 
substitution of the relation (known as ATR scanline) , 

k<^pp — q — ksuvO — —^/e^suvO , (113) 
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into the SPP dispersion relation (63 1. Rearranging Eq. (113), the energy of the 
photon reads 



qhc 



63 sm f 



(114) 



It becomes possible to excite SPPs when the in-plane momentum of the photon 



(113) coincides with the momentum of the SPPs at the same frequency. Under this 



condition, the value of the reflectance at minimum can be represented as 



R„ 







{Uj£l - 


eg ^ 


(we2Ki)^ 




+ k. 




£0 J- / 2 


(^£2^1)2 



(115) 



Notice that the expression (115) must be evaluated at uj — uj 




In Fig. [13] the reflectance and the absorbance are plotted as functions of the 
incident angle 6 and the spacer thickness d. We see that there are some optimal 
values of 9 and d for which the reflectance is almost zero. They correspond to an 
efficient excitation of SPPs in the graphene sheet. 
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d (|im) 



Fig. 13. Absorbance and reflectance spectra. Left: reflectance |rp (solid line), and absorbance 
1 — |rp (dashed line), as function of the angle of incidence 9/Oa- Right: reflectance |rp (solid line), 
and absorbance 1 — |rp (dashed line), as function of the spacer thickness, d. In both panels the 
photon energy m fiw = i.l meV and other parameters are as in Fig. |11| 




Fig. 14. Dependence of the resonant frequency ojniin (a^c), and reflectance at resonance iJmin (d- 
f) upon angle of incidence 9, obtained from numerical calculation [red solid line], approximations 
| |119| [blue dashed-dotted line] and \122\ [green dashed line]. The parameters are: d = 2.51 ^m 
(a,d), d = 7.53 fiva (b,e), and d = 12.55 /^m (c,f). Other parameters as in Fig. |11[ 



6.4.2. The limit K2d <C 1 



Putting = in Eq. (110), wc find 



E. 



(r) 



KiCJEa + ikzA^^^ 



ik,A 



(+) 



g - k^cr'Ki/eo 



(116) 
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where ^ is defined as 



cr 



(117) 



The above equation is the total internal reflection (91 1, modified by the presence of 
graphene. The smallest value of E'i''^ / Ex'' is obtained for Im^ — 0, that is, for 



a" 



0, 



(118) 



which coincides with the dispersion relation of SPP waves with the wavenumber ki 
in a medium of dielectric permittivity ei/2. We recall, however, that the fulfilment 



of the condition ( 118 ) by itself does not imply an efficient excitation of SPPs because 



this type of electromagnetic waves cannot be excited by direct illumination. We note 
that Eq. (118) fixes a relation between u) and 0^ which are no longer independent 



variables. Solving, for example, Eq. (1181 for sin0, Eq. (116) can be regarded as a 
function of cj only. 



For finite n^d <^ 1 Eq. (112) becomes 



\J £3 sin^ 6* - ei 



AaEpd 



ce2 



£3 sm 



£2 



(119) 



The main consequence of Eq. (119) is the existence of a low-angle cutoff, that is, 

(120) 



SPPs can be excited only for incidence angles higher than a certain Q 

he 



sin^ 9„ 



£2 
£3 



1 



AaEpd 



This expression was obtained from Eq. (1191 by setting uj- 



— 0. Decreasing the 

thickness of the gap between the prism and graphene, d, leads to an increase of 
©min- When the right-hand side of Eq. (120) becomes larger than unity, SPPs no 



longer can be excited in the full range of incidence angles 6. 



14 



The validity of approximation ( 119 ) is demonstrated explicitly in Fig. 

=a 1.236'c. At the same time. 



a). For 



small d, the approximation (119) is valid for 6* > 



for 9 > 2.159c one obtains a good agreement between the approximation ( 119 ) [blue 
dashed-dotted line] and exact results [red solid line] both for Wmin [Fig. 14 'a)] and 
for -Rmin [Fig. [Wjjd)]. Finally, increasing d limits the validity of the approximation 



(119) to small 9 only [see Figs. Mb), Mc), pTe), and 14 [f)]. 



6.4.3. The limit K2C? ^ 1 



In the opposite limit of K2d oo, we obtain for Eq. (110) 

^263 - fcze2 



1^263 -I- ^262 ' 



(121) 
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which is equivalent to Eq. (91 ), as expected, because graphene is far away from the 



2-3 interface. Solution of Ec^. (112 1 by means of perturbation theory yields: 

-1 

AaEp El £2/3 

^inin 



£3 sm y — £i 



£3 sm 



£2 



where 



Wo 



— tanh 



£3 sm y - £2 



(2 



£3 sm 



£3 sm 



£2 



(122) 



(123) 



Comparison of panels (a-c) in Fig. 14 demonstrates that increasing d leads to a 



better correspondence between the approximation ( 122 ) [green dashed lines] and 
the exact result [red solid lines]. 



6.4.4. Transfer matrix method 

We will show now how the same problem of interaction of light with a graphene- 
based planar structure can be solved using the transfer matrix method^^. Obviously, 
the result is the same, however, the method is suitable for systems based on graphene 
multi-layers, such as that of Fig. [19) Referring to Fig. flO) the boundary conditions 



(99) and (103) at the 3-2 interface can be written in the matrix form as 





1 


£3 




ZK2 *K2 




'^2,x 




2i£3K2 


.£3 -^2. 




.-£2 £2 . 







^2,x 
^2.2 



(124) 

The determinant of M^^^ is det M^^^ = ifcz£2/(K2£3) • Similarly, the boundary 
conditions at the 2-1 interface can be written as 



^2,x 
^2.2 



1 

2^ 



Xi 
Xi 



1 



1 



(+) 



M 



2^1 



■(+) 

l,x 
(-) 

l,x 



(125) 



where A^*-* has been defined in (109). The determinant of AP^^ is, det AP^^ = 
K2ei/ ('«i£2) • The full transfer matrix is M = M3^2-^^2=^i and its determinant reads: 



det M 



(126) 
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The reflection coefficient is given by 



E. 



-('0 



M12 



(127) 



or explicitly, 



E^"^^ _ K2 cosh(K2d)(xie3 + ikze2^['^^) + sinh(/«2(i)(«;ie3Ai^^ + ik^^e2Xi) 



Ex 



K2 ' 



cosh(K2d)(xie3 - ifc^e2A^+'') + sinh(K2d)(K^e3A^^^ - ik;,e2Xi) 



(+) 



(128) 



where Mij are the elements of the matrix M. Notice, that Eq.( 1271 can be obtained 
from Eqs.(125) and (124), by putting e['^J = which is necessary to guarantee that 
the field is finite at z — >■ 00. After some algebra, Eq. (128) can be shown equivalent 



to Eq. (110) 



6.5. Otto configuration: s— polarized wave 

We shall now discuss the reflection of an s— polarized wave in the geometry depicted 
in Fig. [To] The wave vectors have been defined in Sec. |6.4[ The fields in the three 
regions are given by: 
Region 3 



B, = (i?W,0,i3W)e^^^'-XsW,0,i3W)e^*^'-^ 



(129) 
(130) 



Region 2 



Region 1 



E- 



2,y 



0, i3W)e^^l^'-- + (i?(7j , 0, )e'^l"'-~. 



(-) 

-1 a'--).^ 



B, = {B\-J,Q,B[:]y^''~'-^ 



The boundary conditions at z = 0, d are; 



E. 



(-) 
2,y 



(131) 
(132) 

(133) 
(134) 

(135) 
(136) 



-K2d 



(137) 
(138) 
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We can express the magnetic field components in terms of the electric field ones. 



To that end we use Maxwell's equations (10), and (12), (14) and obtain 







^ z 




rn,x 


UJ " 


5(±) 





These allow to write Eqs. ( 136 ) and ( 138 ) as 



and 



(139) 
(140) 
(141) 
(142) 

(143) 
(144) 



respectively. The problem now reduces to the solution of the following system of 
linear equations: 



Here 



E2,y /Ey 
1 



(145) 



1 

— K.2 



(aP + «;2)e'^2d ) _ ^^y-n2d 
1 



Ki ± ifiQUJcr . 



From the solution of 1 145 it follows that the reflection coefficient is given by 



E., 



('■) 



1^2111,. 



E, 



(146) 

(147) 
(148) 

(149) 



where rji^s = aI ' cosh(K2(i) + K2sinh(K2d), 772, s = As '' sinh(K2(i) + K2 cosh(«;2(^)- 
We note that for c? — J- 00, r has a pole for aI ■* + K2 = 0, which yields the TE-wave 



(-), 



dispersion relation [compare with Eq. (81)]. 

In Fig. [15] the reflectance of a s— polarized wave is plotted against the photon 
energy and the angle of incidence. Clearly, the reflectance is close to unity every- 
where excepting w — > 0. Accordingly, the absorbance, 1 — jr^P, is close to zero. This 
behaviour contrasts with the case of p— polarized waves, where sharp SPP-related 



resonances occur. There are such dips in the reflectance depicted in Fig. 15 (note the 
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0.000 



Fig. 15. Reflectance of a s— polarized wave in the geometry of Fig. |10| The parameters are the 
same of Fig. |11| 

vertical scale) because the conductivity has been modelled by the Drude formula 
(which has a" > 0). Hence, TE plasmon-polaritons in graphene are not captured 
by the present calculation (c.f. Sec.[5|. 

6.6. Otto configuration: arbitrarily polarized waves 



In the general case of an arbitrary linearly polarized wave (see Fig. 16 ), the incoming 
magnetic field is, 

Bi = Bi{— cos (yscos 0, sin (/s, cos Lp sin0)e''^* '' . (150) 

Note that if — Q corresponds to a purely s— polarized wave, whereas if = 1112 
corresponds to a purely p— polarized wave. 

The incoming electric field is obtained from Maxwell's equations as 

Ei — —^ki X B— '^Bilsinif cos 9, cos If, — sin if sin 9)e^^''^ (151) 
k 

and the reflected field can be written as 

Br = {Bs cos e, Bp, Bs sin e)e^^-'^ . (152) 
Again, it follows from Maxwell's equations that the reflected electric field is given 

by 

Er = -^kr X B = '^{-BpCos9,Bs,-BpSind)e'^--^ . (153) 

It becomes clear that the arbitrarily polarized wave can be decomposed into s— 
and p— polarized components, which are reflected independently. We have already 
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Fig. 16. Geometry for the scattering of an electromagnetic wave containing both s— and 
p— polarization components. 



computed the reflection coefficients for s— and p— polarized waves. Thus, we have 
for the s— polarized component: 



cosy'^. (154) 



ijW i?,cos(p' B, 

For Lp = 7r/2, there is no s— polarized reflected wave because the incom ing ER is 
purely p— polarized. Notice that Ey / Ey has been defined in Eq. (149 1. 
For the p— polarized component we have: 

Br, B„ E^J^ , , 

— 7^— — ; ^ = — smos — , (155) 

B,;sin(^' B, ^40' ^ > 



and Ex /Ex has been defined in Eq. ( 110 ). This result implies that for ~ Q there 



is no scattered polarized wave because the incoming ER is purely s— polarized. 

In the regime where the optical conductivity of graphene is dominated by the 
Drude term, the excitation of a s— polarized surface waves in graphene is not possible 
because the imaginary part of the conductivity is positive (refer to Sec.jS]). On the 



other hand, choosing the parameters of the problem appropriately (refer to Sec. 6.4 ) 
it is possible to obtain a strong absorption of the p— polarized component of the 
impinging ER. As a consequence, the setup of Fig. [9] can work as a polarizer. 

If we write B^/Bi = ^^6*''= and Bp/Bi = ApC^^^ , the total reflected power is 
defined as 

Rt = Rs + Rp = \As\^ + \ApW (156) 
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and the angle of ellipticity, tp, is obtained from'Sfil 



2A A 

tan(2V') = ,2 ^ I2 (^os{Sp - 5^) . 



(157) 



The angle vanishes when either Ag or Ap is equal to zero, corresponding to a 
linearly polarized wave, or when 6p — 6s = ±7r/2, corresponding to an elliptically 
polarized wave. 




Fig. 17. Total reflectance Rt and the angle of ellipticity ip as a function of the Fermi energy 
Ep, for different angles of polarization of the incident wave. The reflectance of the p— polarized 
component vanishes for Ep = 0.45 eV. Other parameters as in Fig. |11| The angle of incidence is 
6 = 1.958c and the energy of the incoming light is hw = 5.1 meV. 



In Fig. 17 we represent the total reflectance Rt, given by Eq. (156), and the 



angle of ellipticity ip, as given by Eq. (157), as a function of the Fermi energy 
Ep, for fixed a; and 9 and different angles of polarization, ip, of the incident wave. 
When if — tt/2, the incoming wave is purely p— polarized and the reflectance for 
Ep = 0.45 eV is zero (for the given u and 9). As (p deviates from 7r/2, a finite 
reflectance for Ep = 0.45 eV appears due to the s— polarized component of the 
incoming wave, as can be seen in the left panel of Fig. [17] In the right panel of the 
same figure, the angle of ellipticity (ip) is plotted. We can see that there are two 
points for which ip is zero, they correspond to = at Ep — 0.45 eV and to 
Sp — Ss — ±7r/2 at a sightly different energy. 
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6.7. Otto configuration in a static magnetic field 

We shall now consider the Otto configuration of Fig. [lO] in the presence of a static 
magnetic field Bq = (0,0, Bo)- We take an arbitrarily polarized incident wave, 

(158) 

where we have explicitly separated the p— and s— polarized components. The start- 
ing point to solve the light reflection problem is to write the ER field in regions 
m = 1, 2, 3 in the most general way. We obtain, 

i) Region 3: 

^3 = iEi'\E^\E^^^)e'^'-'' + (i;W, 4'^^ )e^^'-^ (159) 

ii) Region 2: 



E,{r,t)^ (i?W,0,i?W) + (0,i?«,0) 



E2 = {Eir]-Ei^lE^^])e^'^'^^'-^ +{Ei-]^ (160) 
iii) Region 1: 



ik\ 



where the wavevectors ki^r^ 



(±) 



6.4 



(161) 

Similar equations 



have been defined in Sec. 
apply to the time-varying part of the magnetic field B. 

The next step is to write the boundary conditions at each interface. At z = 0, 
the magnetic field is continuous, as well as the tangential components of the electric 



field. The boundary conditions have the same form as Eqs. (99)-(100), and (1351 



(136), and shall not be reproduced again. The boundary conditions at the second 
interface [z — d) must take into account the discontinuity of the magnetic field 
across the graphene layer due to the existence of a surface current Jg [see Eq.([7|)], 



Bi-B, 



IIqJs X n — ^Qa{u!, Bq)Eh (z = d) x n , 



(162) 



with denoting the in-plane components of the electric field with respect to the 
graphene sheet. The magneto-optical conductivity tensor of graphene (j(w, Bq) was 
discussed earlier, in Sec. 3. The boundary conditions relating the electric field am- 



plitudes coincide with those given earher, namely, Eqs. ( 101 ) and ( 1371, whereas for 



the magnetic field amplitudes we obtain the following set of conditions: 



-Mo [crxxE[ J + (TxyE[ J 



B 



(-) 

l,x 



( + )„K2rf 



Mo ( OyxEi 



{-) 



(163) 



(164) 
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As before, we express the boundary conditions using the electric field components 



only. Using Eqs. (103)-(143), we easily obtain 



K2 



I 

£0 

K2 



lujflo 



'^yyEi^y 



(165) 



(166) 



Combining Eqs. (199|), ( |103[ ), ( [135^ , ( |143| ) and (|165[)-(|166|, and defining the following 
vectors. 



(167) 
(168) 



e, ^ (4'^)/£;«,£;W/i?w,i?(7^V4^^), 

a closed system of equations for the ER amplitudes is obtained according to 

where /s > denotes the relative fraction of the incoming electric field stored in 
the s— polarized component, = Eyj/Ex,i, and Sp, $p, S^, $s have been defined 



(169) 



in Eqs. (106), (107), (1461, and (1471, respectively, while other ^-matrices are are 
given by 




(170) 



(171) 



It should be stressed that in Sec. 6.6 the s— and p— polarization components 
were independent, while in the case under consideration (with external magnetic 
field) they are coupled through the matrices S^p and Spg . The optical properties of 
the the system in Otto configuration, in the presence of an external magnetic field 



are fully determined once the system of equations ( 169 ) is solved. For instance, the 
reflectance is obtained according to 

\Ei'^/E^^\^ + cosH\Ei'^^/Ei'^\^ 



R = 



1 + /I COs2 I 



(172) 



and reduces to the familiar form R = {Ex^'' / Ex ^ \ in the particular case of a pure 
p— polarized incident wave (/^ = 0). Note that, by expressing the incident and 
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reflected electric fields in the form similar to Eqs.(151), and (1531 



i?3 ~ (sin iy9 COS ^, cos (p, — sin sin 0)e* 



(— -Ep cos 6', Es 



-Ep sin I 



(173) 

\Es\'^)/\E^\'^, and we have 

/s = tan" 

In Fig.flSlwe plot the reflectance as a function of w for polarization, obtained 



the reflectance assumes the simple form R = (|i?pp 



by solving numerically ( 169 ). We focus on the parameters considered in the absence 



of magnetic fleld (refer to Fig. 13 1. The suppression of the reflectance within speciflc 



frequency intervals, interpreted as the excitation of SPPs for the combined layered 
structure in zero-fleld (see Sec. 6.3), is still seen to occur. However, for i?o > 
the dips are less pronounced and shifted relative to the zero-field case. The dips 
in the reflectance spectrum for non-zero magnetic fleld are related to a special 
type of surface plasmonic waves known as magneto-plasmon polaritons, resulting 
from the hybridization between plasmonic waves and cyclotronic modes. In the 
semi-classical regime, the energy of a magneto-plasmon polariton mode, in the non- 
retarded approximation uj <^ cq P^ l ioo ^ jg given by 



where is the SPP dispersion in the absence of a magnetic fleld for graphene 



cladded by two dielectrics of permittivity £2 and ei [see Eq. (64|], and the cyclotron 



frequency Wj, is given by Eq. ( 33 ) . A similar formula holds for the 2D electron 



gasESmnil Away from the semi-classical regime (i.e., for low Fermi energy or high 
magnetic fields), the quantization of the electronic spectrum must be taken into 
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account. In the quantized regime, the physics of graphene becomes quite exotic, 
with interband transitions between Landau levels in the valence and conduction 
bands of graphene controlling the characteristics of the magneto-plasmon polaritons 
(e.g., their confinement and propagation) and giving rise to quasi-transverse electric 
modes in specific frequency windows. The reader can refer to Ref. for further 
details. 



The formula (174) can be used to estimate the shift of the resonance positions 
due to the plasmon hybridization with cyclotronic modes. It predicts a shift Q{B) — 
n{Q) of roughly 0.56 meV for B = 1 T, and 1.72 meV for B = 2 T, which agrees 
qualitatively with the exact numerical results presented in Fig. fTsl We remark that 



the absolute position of the reflectance dips cannot be obtained from Eq. (174), 
since this result has been derived for a simpler geometry consisting of graphene 
sandwiched between two dielectrics (see also discussion in Sec. 6.3). An interesting 
feature of the hybrid plasmon-cyclotronic modes is the absence of a clear transverse- 
magnetic character (i.e., contrary to SPPs, magneto-plasmon polaritons have non- 
zero magnetic field longitudinal components, and B^^^- Strictly speaking, these 
modes cannot be excited by means of pure p— polarized light, hence explaining why 
the reflectance does not drop all the way to zero at the magneto-plasmon polariton 



frequency (see Fig. 18). For illumination with a TM wave, complete supression of 
the reflectance can only occur when pure SPP (p-polarized) modes are excited (that 
is, at zero external magnetic field). 



6.8. Transfer matrix method for N— layer problem 




We have computed the TM SPP spectrum for graphene double layer in Sec. |4.4[ 
It is instructive to apply the transfer matrix method to a more general problem 
considering an array of N graphene interfaces. It allows for the calculation of the 
reflection coefficient in the ATR configuration whose zeros correspond to the SPP 
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modes (as we will see below, for a propagating wave the poles of the reflection 



coefficient can be related to the SPP modes). Figure 19 depicts the particular case 
of = 2. Thus, the structure under consideration contains + 2 dielectric layers, 
first and iV+2-th layers being semi- infinite with dielectric constants ei and 6^+2- At 
the same time, other layers are characterized by width Dm and dielectric constants 
e™, m = 2,...,A^-|-1. N graphene layers are arranged at interfaces between rn-th 
and m — 1-th layers (to = 2, -I- 1). The final solution can be put in the form: 



(r) 



where the total transfer matrix M 



M 



(+) 



(-) 



M 



m+1- 



(175) 

and matrices M™+i^™ 



M 



can be obtained as generalization of the matrices in Eqs. (124)-(125), namely 

-i [£^+2 kz 

2eAr+2K7V+l [£^-1-2 —kz 

K„ + lD,„ + l Q 

gK^+lA^^ 



lKN+1 



(176) 



M 



m+1- 



1 



2Xr, 



1 

-1 



1 



1 



The determinant of the total transfer matrix is given by 

det(M) = z-^^ 
KieAf+2 

and the reflectance of the structure can be evaluated from the matrix elements 
[Eq. pTl )]. 



,(+) 



(177) 



(178) 



In Fig. 20 we present the reflectance and the absorbance of double-layer graphene 
{N — 2) as functions of both the energy and the incident angle. The most notable 
effect is the shift of the resonant energy from hu w 4 meV to hu 8.5 meV (more 
than 100%!). This is due to the hybridization of the plasmon-polarion bands of the 
two graphene sheets when they come closer to each other. It is possible to control 
the position of the resonance by tuning the distance D2. 



7. SPP excitation by incident wave at a metallic contact 
7.1. Green's function for ER scattering problem 

Here we shall consider the problem of scattering of electromagnetic radiation by a 



very thin metallic stripe (of width L) on top of a graphene sheet (see Fig. 21). It 
has been solved numerically by Satou and Mikhailov i^Sl^ here we shall show how an 
approximate solution can be obtained analytically using Green's function methodic. 
The scattering problem for a plane p— polarized wave leads to the two-dimensional 
Helmholtz equation ^Sil, 

, ,2 



e{z) — By=0, 



(179) 
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0) (meV) 




Fig. 20. Top: Absorbance (dashed line) and reflectance (solid line) versus photon energy of double- 
layer graphene for 9 = 2.16c. Bottom: reflectance as function of both the energy hu) and the angle 
of incidence, d/9c- The parameters are: Ep = 0.45 eV, V = 0.1 meV, £4=14, £3 = 1, £2 = £1 = 5, 
D3 = 2.51 ^lra, D2 = D3, and 6»c = 36.7°. 



where ^ is the Laplace operator in the xz plane, e(z) £2 for ^ < and e(z) = ei 
for z > 0. Considering an incident wave coming from z — —00, the solution of Eq. 



( 179 ) for < can be written as 



By{x, z < 0) = IB, cos(fc,z)e"?" + Bg^at • 



(180) 



The first term in Eq. (180) represents total reflection (as if the whole plane 



2 = were covered with a perfect metal) and the second one (that we call scattered 

field) includes the contribution of all parts of the interface that are not covered by 

the metal (but covered with graphene). As before, here = kcos9, q = ks'mO 
1/2 

and k = uj^^ /c. The scattered field can be presented in terms of Green's function 
G.(r;r-), 



scat 



[e(a;' - L/2) + e(-x' - L/2)] 



{GjX,T-x\z' 
dB,. 



dz 



0)x 
dx' , 



(181) 
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graphene metal graphene 



2 



electric 



Fig. 21. Graphene sheet on a dielectric substrate with a metal contact (stripe of width L) on top 
of it. 



which is given by ' 



G,-(f;r') = ^[i/W(fci?)+i/«(fci?- 



(182) 



with R ^ y^ix-x'Y + {z- z'Y and R' ^ ^{x-x'f + {z + z'Y. Here H^^^ de- 
notes the Hankel function of the first kind and zero order. The Green's function 



( 182 ) is the solution of the equation 



(183) 



and takes into account the effects of a field source located at f* = (x' , z') and of an 



image source at r' — (x', — z'). The Heaviside functions &{x) in Eq. ( 181 1 explicitly 
take into account the fact that the sources of the scattered field are located at 
\x\ > L/2. 



Similar to Eq. (1801, the field in the substrate {z > 0) can be written as 



By{x, z > 0) 



[H^'\[e^/e2]'/^ky/ix - x'y + z^)x 
dB 



[Q{x' - L/2) + <d{~x' - L/2)] 



=0+ 



dx' 



(184) 



7.2. Integral equation 

Matching conditions at z = 0, for \x\ > L/2 read: 



1 dBy 
£2 dz 



z=0- 



1 dBy 
ei dz 



z=0+ 



and [see Eq. (50)] 



By{x,0-) - By{x,0+) 



(j{uj) dBy 



ieQUje2 dz 



(185) 



(186) 



12, 2013 2:29 WSPC/INSTRUCTION FILE UMPByQ 



40 



Substituting Eqs. (pO|, (fTsTl), (pij), and (jlSSl) into (fT86|, we obtain: 

dx' 



1 



262 



2 J|k'|>L/2 



, dBy{x',z) 
g{x -X) 



dz 



z=0- 



cr{uj) dBy{x,z) 



c, ■ (187) 

This is an integral equation for the magnetic field derivative (equal to ifiQUjEx) on 
graphene, with 

gix - x') = e2H^'\k\x - x'\) + e,Hl^\[e^/ e^Y'^k\x - , (188) 

Let us introduce an analytic function /(x), such that it is equal to this derivative 
(with e2Bi = 1) on uncovered parts of the graphene sheet, 

1 dB^ 



£2^* dz 



z=0- 



\x'\>L/2, (189) 
however, not necessarily vanishing on the metal stripe. It obeys the equation, 

gix-x')f{x')dx' , (190) 



SoijJ 2 



\x'\>L/2 



that holds for — oo < a; < oo. Multiplying Eq. (190) by exp(iQa;) and integrating 
we obtain: 



lAnSiq + Q) = ^fiQ) - ilC{Q) 



EqUJ 



e'^"" f{x')dx' , 



(191) 



\x'\>L/2 



f{Q) - / e^'^^fix)dx 



where 



is a Fourier transform of function f{x), and 

POO 



/C(Q) 



K2{Q) Ki{Q) ' 



Ki{Q), K2{Q) were defined in Eq.(44) with replacing q ^ Q. The integral in the last 



term of Eq. (|191|) can be written as 
1 



/ e"^""' f f{Q')e-"^'''' dQ'dx' 

J\x'\>L/2 J-oo 



fiQ) 



1 



1 ^; 



Therefore, Eq. (191) reads as 



47r<5(<z + Q) = i^/(Q)+/C(Q) / CiQ ~ Q')f{Q')dQ' . (192) 



This is a Fredholm integral equation of the second kind ^ with the kernel 

rfn n'^ A^n n'^ sin [(Q - QOL/2] 

^Q - Q ) ^ S[Q - Q ) — — — . 

7r(Q - Q') 



(193) 
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First, let us consider the trivial case of L 
C{Q - Q') = 5{Q - Q') and we have: 



(no metallic stripe). Then 



k{Q) = - 



ATl5{q + Q) 
K{Q) 



.ct(i£) 



and 



/o(x) 



2 exp(iga;) 



(194) 



Substituting (194 1 into Eqs. (181) and (184), we obtain: 



By{x,z> 0) =tj,B,e*(9^^''i(«)^) 



(195) 



with 



]C{q)+i 



2£2 

K2(g) 



/C(g) 
ei 



2 



-i(g)/C(g) + z^' 



(196) 



Here we have taken into account, that ^2(9) — —ikz- We have used the following 
relation for the Fourier transform of the Hankel function: 



dxH^p {k\/ + z^) ex.p{—iQx) — 



2exp(Vfc^ - Q'^x) 



(197) 



Of course, there is no real scattering in this case, just reflection, and the coef- 
ficients Ty and ty could be found in a much simpler way but this simple example 
demonstrates how the method works. Let us notice that, if we formally suppose that 
^1,2(9) are real (respectively, is imaginary), the poles of the reflection and trans- 



mission coefficients (|194|) yield the dispersion relation of surface plasmon-polaritons 

o check that the equation 

/C((7)+^^=0 (198) 



in graphene. Indeed, it is easy to check that the equation 

a{uj) 



eoLo 



is equivalent to Eq. (47) 



7.3. Solution for qL « 1 



We shall now solve Eq. (192) for the case when L is small, QL << 1, where Q is 
any relevant wavenumber along x. Then we put sin [{Q — Q')L/2] ~ (Q — Q')L/2 
and 



C{Q - Q') = 5{Q - Q') - 



L 
2tj: 



(199) 
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By noting that 



f{Q)dq = 27r/(Q) = const , 



Eq. ( 192 1 is solved to give 



/(Q) = /o(g) + /(o) 



IC{Q)L 



/C(Q) 



and 



/(a;) = /o(a;) + /(0) 



L 
2^ 



/C(Q) 



-iQx 



oo IC{Q) + i 



dQ, 



(200) 



where /o (x) is given by Eq. ( 194 1 . The integral in Eq. ( 200 ) converges, since IC{Q) 



Q ^ for large Q, and can be calculated using Jordan's lemma and the contour shown 



in Fig. 22 (for a: > 0). We denote by Qp the pole that lies in the lower half-plane and 



Reg 




Fig. 22. Contour for integral in Eq. |200| for x > 0. For x < it has to be closed in the upper 
half-plane. 



note that it is one of the two roots of Eq. ( 198 1. In other words, Qp and Qp = —Qp 
are two possible values of the SPP wavevector for a given frequency w. Calculation 
of the integral yields 



JCiQ) 



2W-^/C(Q)+^^ 



-iQx 



dQ — i 



ICiQi 



dK(Q) 




dQ 


Q=Qp 



and we obtain: 



Q = 



a{uj) KliQp)KliQp) 

£ow Qp [e2Kl{Qp) + emUQp)] 

fix) = foix) + miQe^'^t- ; 
crju) i4{Qp)i4{Qp) 
£oW Q^ [t2Hi{Qp) + tinliQp)] ' 



-iQ„ X 



(201) 
(202) 
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By requiring self-consistency of Eq. (202), we obtain 

1 2 



/(O) 



and the final result is: 



}C(q)+i 



(203) 



0.05 



10 20 

©(meV) 



03=10 meV 




20 40 

X (jam) 



Fig. 23. SPP excitation by incident ER at metallic stripe. Left: Amplitude A of the second term 
of the function /(x) [Eq. | |203[ )] versus frequency u for different angles of incidence as indicated. 
Right: Spatial profile of the generated SPP wave for two different frequencies indicated on the 
plot. The parameters are: Ep = 0.45 eV, ?ir = 0.1 meV, £2 = 1, ei = 5, and L = 4 fira. 



Equation (203 1 is valid for a; > 0, while for a; < we have — > Qp . The second 
term of this equation represents the SPP wave excited by shining ER radiation at 



the metallic contact. The two waves described by Eq. (203) can be referred to as 



forced and intrinsic oscillations, respectively, of the free 2D electron gas in graphene. 
The metallic stripe is a topological defect imposed on graphene and, as mentioned 
in the Introduction, can be used to launch SPP waves in this two-dimensional 
conductor. The amplitude of the SPP oscillations is proportional to the width of 
the contact ("defect power") and depends on the ER frequency as shown in Fig. 
Fig. 23 Let us remind that the solution ( |203[) is valid in the limit |Qp|L << 1. As 
shown numerically by Satou and Mikhailov SS^ for larger |Qp|L the SPP amplitude 
becomes to oscillate as a function of frequency or contact width. The generated 
polaritons decay with the distance from the contact-stripe antenna because Qp is 



complex (see Fig. 23). The scattered electromagnetic fields above and below the 



interface can be found by substituting (203) into Eqs. (181) and (184), leading to a 
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rather complex dependence on x and z (or, in other words, on the scattering angle 
and distance from the contact). 



8. ER coupling to graphene with periodically modulated 
conductivity 

8.1. Polaritonic crystal 

So far, we have been considering graphene as a perfectly flat membrane with a 
homogeneous optical conductivity. Now we shall relax the latter assumption still 
keeping the former (i.e. graphene's flatness). There are different ways of imple- 
menting a position dependent conductivity in graphene: (i) patterning graphene 
micro-ribbons in an otherwise homogeneous graphene sheet, (ii) inducing an inho- 
mogeneous strain profile in graphene, (iii) producing an inhomogeneous profile of 
adsorbed atoms or molecules on graphene sheet, and (iv) using patterned gates. 
The last three mechanisms produce a position dependent electronic density leading 
to a position dependent conductivity of graphene because of the dependence of the 
Fermi energy on the density. 

We assume a periodically modulated conductivity, a{x + D) = a{x), and in- 
troduce the one-dimensional (ID) reciprocal lattice wave vector, G = 2n/D. The 
system we have in mind is represented in Fig. |24| It can be seen as an ID polaritonic 
crystal rani 




Fig. 24. System with a periodically modulated conductivity period D (side and top views). 



Since the system is periodic, the fields can be written in the form of Fourier- 



Floquet series (in other words, they obey the Bloch theorem). In Fig. 24 we have 
represented an array of graphene micro-ribbons, of conductivity Ug, separated by 
regions of conductivity (Tm, which can be made of either of a dielectric or a metal. 
However, we note that the formalism we develop below applies to any profile of 
modulated conductivity. Thus, Fig. [24| should be understood as a schematic repre- 
sentation of one of the aforementioned possibilities. 



12, 2013 2:29 WSPC/INSTRUCTION FILE UMPByQ 



45 

8.2. Formalism and results 

First, we want to find the band structure of the TM-type surface plasnion-polaritons 
for a system of Fig. [24j If it is composed of periodic patches of graphene, the 
approach presented below fails because it requires an infinite number of reciprocal 
lattice vectors for achieving convergence. (The scattering problem we study later is, 
however, well defined.) 

As before, we assume that the system supports surface waves propagating along 
the graphene sheet and decaying away from it. We consider p-polarized waves writ- 
ten in the form of Fourier-Floquet series, 

BraA^, z)=J2 S„^,||„e^(«+"'^)^e-«'"ii"l^l ; (204) 

n 

EmA^, z)^Y. ^m,.||ne^^^+"'^)^e-'^'"NJ-l ; (205) 

n 

Em.A^,z) = ^£:™,,||„e^(«+"'^)^e-«'"N"l-l , (206) 

n 

where the index m labels the media above (to = 1) and below (m = 2) the z = 
plane, and the sum runs over all integers. Since the system is linear, each Fourier 
component is independent of the others. Then, it follows from Maxwell's equations 
that 

i^7n\\7i = V(g + fiG)2 - w2e„i/c2 , (207) 

^m,y||n = (^1)™-^ —Sm,x\\n, (208) 



fm,2||n = (-1)'"'''^*^— £m,x\\n- (209) 

'^m 1 1 n 

Since the conductivity is periodic, it can be expanded in Fourier series as 

a(x) = ^e^"=^az, (210) 
I 

where 

1 

a, = - / a{x)e"^^'^dx. (211) 
D Jo 

The boundary conditions Ei^x{x,0) — E2^x{x,0) = 0, Bi^y{x,0) — B2.y{x,0) — 
—^oO'{x)Ei^x(xA) read: 

'^^{^i,x\\n — £2,x\\n)e^"''^^ — ; (212) 

n 

- S2,,||n)e"^" = -MO E S^i^i,x||pe^^'+^"'" > (213) 

or, alernatively {I + p = n ^ I = n — p), 

^(^i,y||n - 'S2,y||„)e'"'^^ = ^ 5„_p£i_j.||pe™'^^ . (214) 

n n,p 
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As a result, Eqs. (212)-(213l imply that 



= 



P 



After some algebra, we obtain a non-linear eigenvalue problem, 



^l.xlln H ^l,x\\n "I fn-p^i 

Kl\\n ■ K2||n WEq ^ 



slip 



= 0, 



(215) 
(216) 



(217) 



for the frequency uj. If we consider only the Fourier component p = n in the sum in 



Eq.(217), we obtain 



+ 



£2 



-do = , 



(218) 



which is just the equation for the SPP frequency of the wavevector q + nG [compare 



with Eq. (47l]. The quantity ctq is the average of the conductivity in the primitive 
cell. 



era = 



D 



a{x)dx . 



(219) 



The presence of harmonics with other Fourier components of a{x), ui with ^ 7^ in 



Eq. (2171 gives rise to the band-gap structure of the polaritonic spectrum, equivalent 
to the so called empty lattice approximation of electrons in a periodic potential 
It can be presented either in the extended scheme if we consider the SPP wavevector 
varying from —00 to 00, or in the reduced scheme if we limit the wavevector to the 
first Brillouin zone, — ^ < q < |j. In the latter case, folding the dispersion curve 
into the first Brillouin zone produces upper branches, i.e. the polaritonic crystal 
structure for the SPP frequencies. 

We note that the non-linear eigenvalue problem can be transformed into a linear 
one if we make the non-retarded approximation. As discussed in Sec. |4j we do not 
expect to obtain an accurate solution for the band at wave vectors close to zero. 
In the non-retarded approximation we have Krn\\n ~ k + nG\. Using Eq. (51 1 and 
writing the spatial-dependent conductivity as 



a{x) 



i — s{x) , 



the Fourier transform of <j{x) can be written as 



(220) 



(221) 



where Si is the ^-th Fourier component of s{x). This allows to write Eq. (217) as 



ei + £2 
\q + nG\ 



^l,x\\n 9 / ,^n—p^l,: 



■\\P 



0, 



(222) 
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or as 

iy\q + nG\ 
2eeo 



'5„_pf 1 — w^i^i ,r||„ , (223) 



which has the standard form of a Unear eigenvalue problem. However, the band 
structure derived from such procedure is not quantitatively accurate. 

For obtaining the spectrum in the vicinity oi q = G/2 = tt/D, we need to take 
into account field harmonics £i,x\\n with n = —1,0 only (which correspond to the 



lowest gap in the spectrum). In this case, Eq. (2231 reduces to two equations only: 

■^r^ {So£l^x\\-l + S-l£i^x\\o) = , (224) 

7^ (^if + i.x||o) — ^'^£i,x\\Q , (225) 

4eeo 

which form a 2 x 2 eigenvalue problem. The eigenvalues are 
2 ( ^ I 1 ^ ^ — \ aEpcG 



4ee( 



(5o ± = (So ± VS^^ . (226) 



Equation ( 226 1 describes the form of the spectrum in the vicinity of the gap at 
the Brillouin zone edge. As expected, the value of the gap depends on the Fourier 
components cr±i of the conductivity. The approximated description of the spectrum 
that we have presented has to be checked against a full numerical calculation. We 
have verified that for the case of an array of micro-ribbons this simplified description 
fails because a large number of reciprocal lattice vectors are necessary to describe 
the spectrum accurately. In fact, this particular case is poorly convergent. This is a 
consequence of the non-continuous nature of graphene in the micro-ribbon structure. 

For further progress we need a model for a[x). We assume a conductivity profile 
of the form 

a{x) = aDs{x) = cr£,[l - /i cos(27ra;/D)] , (227) 
where ud is given by Eq. ( [23^ . The Fourier transform of s{x) reads: 

50 = 1 , (228) 

51 = -^(5iA + 5i,-i). (229) 



For this form of s{x) the eigenvalue problem (2231 has a fast convergence. 

In the top plot of Fig. [25]we present the band structure of the polaritonic crystal 
produced by the conductivity modulation, computed using the non-retarded approx- 
imation. The dashed black lines represent the folding of the bare SPP dispersion 
curve of homogeneous graphene into the first Brillouin zone. The solid lines repre- 
sent the spectrum when the conductivity is modulated. Clearly, there is a large gap 
between the first and the second bands and smaller gaps between the other upper 
bands. The blue dashed straight lines represent the light cone, hckj ^ft{. The states 
of the upper bands located within the cone can be excited by shining light on the 
flat graphene sheet without the aid of a prism. 
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q (1/|^m) 



Fig. 25. Band structure of the polaritonic crystal obtained from the solution of the linear eigen- 
value problem | |223| , i.e. in the retarded approximation (upper panel) and by directly solving Eq. 
| |217| (lower panel). Only the first six bands are shown. The dashed curves in the upper panel 
correspond to the "empty lattice" approximation. Notice the anti-crossing of the bands n = ±1 
in the vicinity of g = (near 15 mev). The parameters are Ep = 0.45 eV, D = 10 fim, ei = 3, 
£2 = 4, and h = 0.3. 



As it can be observed in Fig. [25j the band structure for is inaccurate k Ri 0, 
the dispersion curve for the lowest band hes inside the hght cone and the curves 
representing the upper branches do not show the expected anti-crossing. These are 
artefacts produced by the non-retarded approximation used. The band structure 
calculated taking into account the retardation effect ^ is free from these artefacts 
(see the lower plot of Fig. 25 ). 
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8.3. Scattering by a polaritonic crystal 
8.3.1. Formalism 

Having studied the emergence of the polaritonic band structure in a system with 
periodicaUy modulated conductivity, we now want to discuss the scattering of light 
from a structure of the type represented in Fig.[24j We consider TM-polarized waves 
in the form 

B^,y{x, z) = <5™,24')e'«"e*'^^ + S™,y||«e^^'+"''^"e-'='"N J-l , (230) 

n 

z) = SmaE^^e^'^-e^'^'' + ^ £„,,,||„e^(«+"'^)^e-«"ii"l^l , (231) 

n 
n 

where q = ksinO, — kcosO, k — ^ye^{uj/c). The relations between the amplitudes 
of the reflected and transmitted fields, Bm,y\\m £m,x\\m and £m,z\\n are described by 



Eqs.(208) and (^209|, while those between the amplitudes of the incoming field are 



the same as Eqs. (|83| and (84), with £3 being replaced by £2 



Notice that in Eqs. (230 ), (231 1 and (232 1 the ratio {q+nG)/k can be interpreted 



as the sinus of the scattering angle of the Fourier mode n, that is, 

^i{q+nG)x _ gifcsine2||„x (233) 

leading to 

02\\n = arcsin[(g + nG)/k] . (234) 

As a result, the usual condition for the Bragg scattering reads as 

A:sin02||n = ksinO + n2'ri / D < k ; 

L»(sin e2\\n - sin 0) = nX , (235) 

where A = 27r/fc is the wavelength of light in the top dielectric with £2- In fact, {q + 
nG)^ — [^121171]^ = k'^ for both propagating and evanescent waves, but for diffraction 
orders (n) corresponding to propagating waves 

K2||n = -«fccos6l2||„. (236) 

In a similar manner it is possible to introduce the scattering angles for transmitted 
waves, 6*1 II „ 

g + nG = (a;/c)£y^ sin6'i||„ , (237) 

f^iWn = -ii^/c)eY^ cos6'i||„. (238) 
The boundary conditions imply: 

fi,x||o =f2,.||o + -B«, (239) 

■^^i,:.!!,! = '^^2,^||n A n^O, (240) 
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and 



P 

Bl^y\\n - B2,y\\n = -fJ'o'^^n-p£l,x\\p A 71 ^Q. 



(241) 
(242) 



Using the relations between the fields, Eqs. (83) and (208), the set of boundary 
conditions reduces to 



£2 



l^l\\0 K2||0 



^l,x\\0 H 2^ <^-p<>l,x\\p — 'l^^x 



UJ€o 

,^l.a;||riH f^n-pSl.xllp — A Tl^O. 

We recall that 

Km\\n = V (k sin 9 + nGY - uj'^e,yjc^ 
for non-negative arguments of the square root, otherwise Hrn\\n is written as 



Km||n = -«V^^em/c2 " (fc siu 6* + nG)^ 



(243) 
(244) 

(245) 
(246) 



The last equation imposes ki||o — —ik^. The choice for the sign of the square root is 
dictated by physical reasons: we must have reflected waves for z < and transmitted 
waves for z > (so called Rayleigh conditions). 



8.3.2. Reflectance and transmittance efficiencies 



Let us compute the fraction of the ER energy carried by the different diffracted 
orders. The time average of the Poynting vector, (S) = ^{E x H*)/2, is given by 
(recaU Sec. 



(5) = ^n{~xE,B;+yE,B;). 
2^0 " " 



Thus, for diffraction order n we have 

5W 



2/ioc2fc^ 



'2,z||r 



We2 



^1, 



:||r 



2/Xoc2|k2||„ 



\E. 



2Moc2|Ki||nl 



Ei^x\\ 



2,a:||n| 
2 



(247) 

(248) 
(249) 
(250) 



for the incident, reflected, and transmitted power per unit area (along the z direc- 
tion), respectively. Here we have assumed that k„i||„ is pure imaginary; if it is real, 
then the corresponding order n carries no energy due to evanescent character of the 
corresponding wave. The negative sign of the reflected power corresponds to the 
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wave propagation in the negative direction along the axis z. Finally, the reflectance 
and transmittance efficiencies are 



7^„ = - 





n kz 




l«2||«| 


\.z\\n 






£2 h^llln 



I '^x\ \n I 
I |n I 



(251) 
(252) 



respectively. The diffuse reflectance and diffuse transmittance amplitudes of the 
order n entering these relations are defined as 

£9. 



c||r 



tx\\n 



■'2,x\\n 



J^x 



£l,x 



11^ 



The specular reflectance and transmittance {n — mode) are given by 



T^a = \rx\\o\' 



and 



To 



^x\\0\ 



|2 €1 



cosO 



£2 



ei/e2 



(253) 
(254) 

(255) 
(256) 



The last expression is valid for ei > 62 or for 9 < 9c = arcsin(\/ei/e2), with 9c 
denoting the critical angle for total reflection at the interface. 

In general, for n ^ 0, using the scattering angles for diffuse reflection and trans- 



mision modes, Eqs. (236), (238), the diffuse reflectance efficiency is 

cos 9 , 1 2 



Tin 



cos 9' 



2||r 



' a;| |n I 



and the diffuse transmittance efficiency is given by 

[ei cos 9 I. 1 2 



Tn = 



£2 COSt/i||„ 



Ik 



(257) 



(258) 



It should be noticed that Eqs. (257) and (258) arc valid for positive arguments of 
the square root only (i.e. for propagating waves) and for real dielectric permitivitties 
ei and €2- The absorbance is defined as A — I ~ J^ni'^n +Tn), where the sum over 
n is restricted to the diffraction orders corresponding to propagating waves only. 



8.3.3. Two special cases 

Let us consider a special case of vanishing conductivity ("graphene is absent"). 
Then Eqs. (243) and (244) give £i^x\\n = (n 7^ 0) (no diffuse reflectance and 
transmittance in this case) and 



£l,x\\Q 



E. 



2 1 



1 - e2 sin^ 9/ei 



(259) 
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The reflection coefficient is given by 



i.^llo ^ y/l - 62 sin^ 6 lei - ^ ei/ea cos d 
Y^l - £2 sin^ 61/ei + ^/t^ft2 cos d 



(260) 



which reproduces the weh known resuh from elementary optics. 

Another particular case is obtained when dg = D. In this case, only the to = 
Fourier component of (t(x) survives. Thus, we obtain: 



Si. 



x\\0 



1 



ei/e2 



p(0 

-Cjx 



COS 6 



ei/e2 - sin^ 



cos 9 



(261) 



Again, using Eq. ( 256 ) and putting ei = £2 = 1 we obtain the well known result for 
the transmittance of free-standing graphene (for a TM wave). 



2 + ajj cosO/eoc 



(262) 



8.3.4. Scattering by a grid of graphene micro-ribbons 

After checking two special cases of the general expressions, we now want to address 



an example where the full set of Eqs. ( 255 ) - ( 258 ) has to be used. This particular 



problem was considered in Ref. For an array of graphene ribbons, the Fourier 
coefficients in Eq. (221) are given by Sq = dg/D and 

lirdn 



Si 



1 



D 



^a_^-il7rdg/D 



I 7^ 0. 



(263) 



In Fig. 26 we represent the absorbance and the transmittance of a grid of 
graphene ribbons for different values of the broadening F. This calculation required 
600 reciprocal lattice vectors. The same quantities for uniform (pristine) graphene 
are represented by dashed curves. For small values of F, both the absorbance and 
the transmittance spectra present a set of resonances, with the most prominent one 
at low energies. On the other hand, for large F, the weaker resonances are washed 
away and the most prominent one becomes less pronounced. Compared to the case 
of pristine graphene, the absorption is suppressed at low energies whereas the trans- 
mittance is increased. The energy of the largest resonance is approximately given 
by hflpi^n/dg), suggesting a coupling of the impinging ER to the SPPs modes of 
graphene. Unfortunately, the analysis is not so simple, because each ribbon has a 
finite size and the electron confinement has to be taken into account. The behaviour 
of the transmittance at low energies should also be noted. In the case of grid it tends 
to unity whereas in the pristine graphene case it attains its lowest values. The grid 
transmits more at low energies although the size of the gaps between the ribbons is 
much smaller than the wavelength of the incoming radiation. 
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(1=0.45 eV, e =3. e,=4, D=8 fim, dg=D/2 




0'- *^ , \ , \ , — I 0' — • — ' — • — ' — • — ' — • — ' 

10 20 30 40 10 20 30 40 
0) (meV) CO (meV) 

Fig. 26. Absorbance (left) and transmittance (right) of a grid of graphene ribbons for different 
values of the broadening F. The dashed lines are the absorbance/transmittance curves for pristine 
graphene for the same F values and chemical potential Ep. The parameters are: Ep = 0.45 eV, 
D = 8 fira, dg = D/2, ei = 4, and e2 = 3. 




Fig. 27. Absorbance of a grid of graphene ribbons: (a) dependence on the lattice parameter D; 
(b) dependence on the angle of incidence; (c) dependence on the electronic density; (d) dependence 
on the width of the graphene ribbon (keeping the width of the unit cell constant). The parameters 
for the reference curve (solid blue) are: Ep = 0.23 eV, D = 8 /im, ei = 5, £2 = 3, and F = 2.6 meV. 
In panel (c) rie is given in units of 10^^ cm"'^. 

In Fig. [27] we study the effect of different parameters of the problem on the ab- 
sorbance curve. In panel (a) we find the variation of absorbance with the parameter 
D, keeping dg — D/2. The shift of the absorbance maximum scales as 




(264) 
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In panel (b), the dependence of the absorbance on the angle of incidence, 9, is 
given. No appreciable effect is seen here, except close to grazing incidence. We will 
see below that in the case of a continuous sheet with a modulated conductivity the 
situation is quite different. In panel (c), the dependence on the electronic density is 
shown (the values on the plot are given in units of 10^^ cm~^). The red shift scales 
with Ue as 



(265) 



Finally, in panel (d) we present the dependence of the absorbance on dg, keeping D 
constant and equal to 8 ^m. The blue shift follows the scaling relation 



Wmax oc 



(266) 



In this latter panel the absorbance of an infinite graphene sheet is also plotted. 
Clearly, there is an enhancement of absorption due to SPP around certain specific 
frequencies, where the absorbance is higher for the grid of ribbons than for an 



infinite graphene sheet. Relations (264), (265), and (266) have the same functional 
form of dependence on dg and Ue as the SPP frequency in the continuous system for 
the wavenumber q — /dg. Considering the case of Fig. 28 we predict a maximum 
at the frequency 



Mir 



-Ephcq = \l— — ;:0.23 X 0.2- 



.1 meV, 



ei + £2 ' ' V 137 X 
a value close to the position of the first maximum in the absorbance spectrum. 



(267) 




0.003 0.006 0,009 0,012 O-OI.-i 0,018 0,021 0,024 0,027 0,03 
to(eV) 



Fig. 28. Contribution of the first three harmonics to the absorbance curve. The reference curve 
(top panel) has F = 0.6 meV. Other parameters as in Fig. |27| In the bottom panel, the quantity 

l^^i,2!|ln/-Ei;''P is depicted. 
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In order to understand the origin of the second maximum (and the third, as 



well) seen in the absorbance spectrum, we plotted in Fig. 28 the squared absolute 
value of the amplitude ^i^xWnl e'^^ for n — ±1,±2 and ±3. Clearly, the spectral 
weight associated with different harmonics is centered at the maxima observed in 
the absorbance spectrum. This representation gives a qualitative understanding of 
how the spectral weight associated with the different SPP modes in the infinite 
sheet is redistributed to form the absorbance spectrum in the periodic system. 
Different harmonics contribute differently to these maxima, for example, none of 
the depicted harmonics contribute to the maximum observed above 24 meV, only 
those with Inl > 4 contribute to it. 



8.3.5. Scattering from graphene sheet with cosine-modulated conductivity 

We next consider an example where the graphene sheet is continuous (as opposed 
to the grid of ribbons) and possesses a periodically modulated optical conductivity. 
As in Sec. 18. 21 we assume a conductivity profile of the form (227), with the Fourier 



harmonics given by Eqs. (228) and (229). We shall see that, for angles of incidence 



different from zero, the absorbance peak will split into two, associated with the 
lattice vectors G = ±2tt/D. We give below a derivation for the energy of the peaks 
as function of the angle of incidence. In the non-retarded approximation, the energy 
of the plasmon-polaritons is given by Eq. ( 64 ), with e = (ei + e2)/2. For an ER wave 
incoming at an angle of incidence 9, the wavevector of the excited SPP is given by 



(268) 



q = |A:sin6' + nG\ — \ — - — sin6' — mG\ 



From Eq. ( 63 ) for the SPP dispersion relation it follows that 



We assume that G > There are two possibilities: 

(1) n = -1 



(269) 



(2) n = 1 



(Sw)^ ei + £2 

Aa EpJt^ 



+ hio sin Q — Ghc/ = ; 



(270) 



(huj)'^ ei + 62 



Introducing the parameters 



hu! sin 6 — Ghc/y/e2 ~ . 



1 £1+62 



AaEp JT2 



(271) 



(272) 
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and 



the solution of the two equations for uj{6) reads as 



b = , (273) 



huj = ±— + — Vsin^ d + Aab, (274) 
2a 2a 

for n = 1 and n = — 1, respectively. Thus, we expect to observe a peak splitting for 
^ > at the left {n — —1) and at the right (n = 1) of the single peak at 9 = 0. We 
shall see that this is indeed the case. 



Ep=0.45 eV; r=0.6 meV; e =3; £,=4; k=0.3; D=1 [xm 




10 12 14 16 18 20 

CO (meV) 



Fig. 29. Spectra of a graphene sheet with a cosine-modulated conductivity. Reflectance (top) 
and absorbance (bottom) versus frequency for different angles of incidence. The parameters are: 
Ep = 0.452 eV, D = 10 fim, ei = 4, £2 = 3, T = 0.6 meV, and h = 0.3. The Brewster angle for 
these dielectrics is 49.1°. 



The linear system defined by Eqs. (2431 and (244 1 is solved numerically and the 
sums over n are cut off at n = —N,...,0,...,N. The numerical solution rapidly 
converges with N. Results for the specular reflectance, TZo, and for the absorbance, 
A = 1 — TZq — To, are given in Fig. |29| For a modulated conductivity the momentum 
of the SPPs is conserved up to a reciprocal lattice vector nG, with n = ±1, ±2, . . ., 
that is, 



Qp ^ \k sine + nG\. (275) 

In this case, even for normal incidence, it is possible to excite SPPs. We stress that in 
the present case the Bragg scattering mechanism expressed in Eq. (275), allowing 
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to overcome the momentum mismatch between the propagating wave and SPPs, 
is induced entirely by the conductivity modulation and it is not a consequence 
of a external grating. The excitation of SPPs at normal incidence is illustrated 
in Fig. 



29 



The dashed black curve represents the behaviour of the system for a 
homogeneous conductivity and impinging ER at normal incidence; clearly the curve 
is featureless. For the inhomogeneous case, a large enhancement of the absorbance 
is seen around the energy given by Eq. (64) with qp = l-KjD (note that using this 
equation implies the extended band scheme). The position of the peak does not 



coincide exactly with the number given by Eq. (64) because this equation is not 



sensitive to the details of the band structure. From Fig. [25] we can observe that the 
bands for n = ±1 at the zone center have an energy of about 15 meV (the energy 
for which the reflectance curve has a maximum for = 0). As the angle of incidence 
approaches the Brewster angle, 0b, for the two dielectrics, the reflectance decreases 
substantially. Note that the Brewster angle of the system is not given exactly by 
the usual formula, 8^ = arctan \/ei/t2, because of the presence of graphene. For 
incidence angles above Qb the reflectance develops two dips and can be larger than 



it would be for 9 — Q (see curve for 9 — 80° in Fig. 29). 

When the incoming beam deviates from normal incidence (i. e. ^ ^ 0), there 
is a peak splitting both in the reflectance and in the absorbance curves, as pre- 
dicted above. We would like to understand in qualitative terms the behaviour of the 



splitting as function of 9. Using Eq. (275) in Eq. (64 1 yields 



-EFhc\ksm9 + nG\ 



(276) 



Clearly, when n > and 9 increases, the frequency shifts toward higher energies. 
On the other hand, when n < 0, the energy decreases as 9 increases. This behavior 



can be understood from the analysis of Fig. 25 For 9 — the light line, Eq. (113) 



is vertical and touches the second band at the center of the Brillouin zone (q = 0) 
where the branches associated with n = 1 and n = —1 almost touch each other 



(they do touch within the non-retarded approximation, see Fig. 25). As 9 grows, 
the slope of light line decreases and the branches with both n = 1 and n 
away. 



-1 split 



In Fig. 30 we plot the two branches of Eq. (274) and compare them with the 
positions of the absorbance peaks obtained from Fig. [29| The agreement is only 



qualitative because, first, Eq. (274) is derived from a kinematic argument and. 



therefore, misses the dependence on h (and, eventually, some symmetries a partic- 
ular problem may have) and, secondly, we used the non-retarded approximation. 
However, for small h the agreement is quite good. Indeed, if we shift the solid curves 
by adding a constant, they would fit the points (solid squares) obtained 

for different 



in Fig. 30 



from Fig. 



29 



30 



we plot |r^||±i| 



as a 



In the bottom panel of Fig. 
function of the energy. The energies of the peaks of |7'a;||±iP are the same as those 
of the transmittance minima (absorbance maxima). This shows that the polaritons 
of the n — ±1 branches with the wavevector gp « are responsible for the features 
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Fig. 30. Position of the absorbance peak as function of the angle of incidence. Top: dependence of 
the absorbance peak energy upon the angle. The lines with squares are obtained from the bottom 
panel of Fig. |29| the solid lines are the two branches of Eq. | |274[ l. The bottom panel shows the 
squared absolute values of the amplitudes of the modes associated with the SPPs of momentum 
itG. The parameters are the same of Fig. |29| 



in the reflectance, transmittance and absorbance spectra. 



9. Scattering of ER from corrugated graphene 
9.1. Setting the problem and definitions 

In Secs.|6j[7]and[8]we have seen that resonant couphng between ER and SPPs can be 
achieved either using the ATR scheme (Otto configuration), or due to a topological 
defect on graphene or its modulated conductivity. Another possibility is the use of 
a dielectric grating ^^^31 jjj -vvhat follows we discuss the coupling of ER to SPPs in 
graphene-based gratings. 

We want to solve the scattering problem of light impinging on a diffraction 
grating, as schematically represented in Fig. [3l] The region 1 is such that z > 
max[a(a:)], whereas the region 1~ is such that a{x) < z < max[a(a;)]. A similar 
definition applies to the regions 2 and 2+ when z < min[a(a;)] and min[a(a;)] < z < 
a{x), respectively. In general, the calculation of the fields in the regions 2+ and 1^, 
i.e. min[a(a;)] < z < max[a(a;)] is a challenging problem. 

It will be assumed that in both regions 1 and 2 the dielectric functions, ei and 
£2, are constant. For a TM wave we have Maxwell's equations (11), (13) and (15), 
while the electromagnetic fields in regions 1 and 2 are the same as represented by 
Eqs. ( 230 ), ( 231 1 and ( 232 1. The Rayleigh-Fano approximation (also called Rayleigh 
hypothesis) assumes that Rayleigh expansions (see below) are valid in the regions 1~ 
and 2+, which is not true in general El (for example, the method fails for rectangular 
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gratings). However, in the case of smooth grating profiles it can be proved that the 
solution of Maxwell's equations exists and is uniqueE^MI^H and can be approximated 
by linear combinations of reflected and transmitted waves in the regions 1~ and 2"*" 
(see Ref. for details). 

Since we are dealing with a corrugated surface, we need to defined tangent and 
normal vectors to the curve z = a{x) (see Fig. 31 ). The unit vector tangent to the 
curve z = a{x) at point x is given by 



t = 



1 



The normal unit vector to the curve is 



da 
dx 



\t\ = l. 



1 



K dx f 



da 
dx 



X + z 



(277) 



(278) 



and clearly ii ■ i — 0. Since the current density vector is tangent to the graphene 
sheet, we must have Jg — uEt — oEft, where Et is given by 



Et=E-t 



1 



Ecr 



da 
dx 



-E, 



(279) 



One of the boundary conditions ([t]) requires the determination of the vector product 
Jc X fi: 



X h = — 



„ da „ , 
Ex + —E, I z . 
dx 



(280) 



and the respective boundary condition Bi,y{x, z)y — B2.yy — ij-quJ x h reads 



Bi^y{x, z) - B2. 



V dx J 



Err 



da 
dx 



-E, 



The second boundary condition is Ei t = E2.t, which can be written as 

Ei^^{x, z) + '^Ei^z{x, z) = E2^x{x, z) + '^E2,z{x, z) . 
These boundary conditions arc to be applied to a general profile a{x). 



(281) 



(282) 
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9.2. Rayleigh method by Toigo et al. 



We give a representation of Rayleigh equations based on Ref. ^. Using Eqs. (230) 



(232), the boundary conditions, Eqs. (281)-(282), can be written in terms of the 



x-component of the electric field as 



E 



l + i 



da q + nG 



dx 



£^^^^^^e'^ii+'^c;)x-^,H„a(x)] . (283) 



E£ ^i[{q+nG)x-K.m„a{x)] 

n 



l+i 



1 - i 



da q + nG 

dx Kl||„ 

da q + nG 



+ 



da q 



dx K2||n 

^{qx+k^a{x)\ ^ (284) 



Eqs. (283), (284) were obtained using Eqs. (83)-(84) (replacing 63 by £2) and (208) 



(209). The boundary conditions are applied at z = a{x). Multiplying Eq. (283) by 



e ''•■(<i^vG)x integrating over dx/ D, we obtain 



l^l\\r 



H2\\n 



2||P- 



k ^ 



. cr / q + nG 

eo^ „ V '«l||n 



(285) 



where 



m\ \p—n 



N — 



D 



fl^Q-iip-n)Gx ^{-1)"' K„^Hr^a{x) 



1 

D 
1 

D 
1 

1 f° , da e-»(P-")G^ 







D ^ — i{p—n)Gx 
. /t , (da-] 



dx- 



dx 



-Ki||„a(a;) 



.. (P-»)G 

'*l||n 



(286) 
(287) 
(288) 



(289) 
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After some manipulations we arrive at the following equation: 
^( ei a iq + nG){q+pG)-uj\i/c^ \ 

n V 111" ° 111" / 

-e2 E ^A^2|b-„ + *^i^«M(^) , (290) 

„ '*2||ri 

In a similar manner, multiplying Eq. (|284[) e^*'^^^''-'^ and integrating over 



dx/D, we obtain: 



E 

n 

E 



(g + nC;)(g+pG') -c^'ei/c' 
{q + nG){q + pG)-uj^e2/c^ 



£l,x\\nMillp^n 



^2.,x\\nM2\\p-n ' 



(291) 



Solution of the boundary problem (2901 and (291) will allow for the calculation of 



the reflectance, transmittance, and absorbance spectra. 



9.3. Three limiting cases 



First, we will show that, if a{x) = and Ex'^ = 0, one obtains from Eqs. (290) and 



(291) the usual SPP spectrum folded into the first Brillouin zone. When a{x) = 0, 



we have Mm||p_„ — 5p_n and Eqs. (290) and (291) reduce to 



r * ^1.2;||p 



2.a:||p 
K2||p 



Solving for E2^x\\p, we obtain 



£2 



0. 



(293) 



'^IIIp '«2||p £01^ 

Since the SPP amplitude decays away from the graphene sheet, we consider k„||„ in 



the form (207) with q e [-tt/D, tt/D]. With this choice we recover Eq. (47) folded 
into the first Brillouin zone. 

When Ex is finite, we have a scattering problem. In this case Mp = Spfl and 



Eqs. (290) and (284) reduce to 



r I '^l.ajllO 



^2.x\\0 , . 
'€2 + I 



«2||0 

^i,2;||o = E; 



2,x\\0 ■ 



(294) 



In the limit c = we obtain Eq. ( 260 ) , the well-known refiectance amplitude from 



elementary optics. In the case a ^ and ei = 62 we obtain Eq. (262), which gives 
the transmittance amplitude for graphene. 
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9.4. A non-trivial example I: sine profile 

The last two examples are trivial for they refer to the limit of zero curvature. We now 
consider a non-trivial case where graphene has a well defined periodic corrugation. 
We assume a profile of the form 

a(x) = /isin(27rx/i:>) , (295) 

from which follows da/dx = {2TTh/ D) cos{2ttx/ D). For this choice of profile the 
Rayleigh hypothesis is exact. We have: 

M^*) = Jj,{k,h) , (296) 
M„||p_„ = Jp-„[*(-l)"+i«™||„/i] = ^''-''Ip-n[{-^r^'^,n\\nh] , (297) 

where Jp{z) and Ip{z) are the usual and modified Bessel functions of order p, re- 
spectively. 



D=10 ^m; e=0; r=2.6 meV; |i=0.45 eV; e,=1 ; e,=1 1 




0) (meV) 



Fig. 32. Dependence of the reflectance, transmittance and absorbance on frequency for different 
depth of the grooves for sine profile gratings. We compare the Rayleigh method with the results 
from the extinction theorem. In the legend, the label R stands for the Rayleigh approximation 
whereas the label ET stands for the extinction theorem result. The parameters are Ep = 0.45 eV, 
D = 10 lira, ei = 1, ea = H (silicon), T = 2.6 meV, and 6» = 0. 



In Fig. 32 we represent the efficiencies TZq and To (top and central panels, respec- 
tively), and the absorbance, A = 1 — TZq — To, as functions of the incoming photon 
energy, for different values of the ratio h/D. When h/D 1, we recover the proper- 
ties of a flat graphene sheet. We compare results from the Rayleigh approximation 
described above with those produced by the extinction theorem method (dashed 
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curves), an exact integral-equation approach to the scattering problem EEH We see 
that the agreement is excellent. In Fig.[33]thc quality of the agreement between the 
two methods is more evident. In this figure we compare the dependence of the spec- 
ular reflectance TZq on frequency, for several angles of incidence 6 and h/D = 0.1. 
The resonance seen in Figs. [32] and [33] above 10 meV is due to the excitation of a 
surface-plasmon-polariton of energy 

hilp = W Epch— « 11 meV . (298) 

V £i + ^2 D 




Fig. 33. Left: Dependence of the reflectance on the angle of incidence for a sine profile grating. We 
compare the Rayleigh method (solid circles) with the results from the extinction theorem approach 
(solid lines). The arrow indicates the direction of growth of the angle 6. Right: absorbance for 
different values of the broadening F, at normal incidence. Other parameters are as in Fig. |32[ 

In the right panel of Fig. J33]we present the absorbance, A, for different values of 
r. As r decreases, the resonance in the absorbance becomes more prominent. Also, 
the coupling of the ER to the SPP of wave number iir/D (extended band scheme) 
becomes evident as a smaller resonance. 



9.5. A non-trivial example II: sawtooth profile 



In the case of a sawtooth profile the function a{x) reads (see Fig. 34 1 



a{x) = 



2h, 
- ) 'D- 



2h^ 



h 
2 ' 



-D/2 < X < 0, 
< X < DI2. 



(299) 
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substrate 



Fig. 34. Periodic patterned substrate with a sawtooth profile. In regions with the label 1 the 
strain is expected to be higher than in those with label 2. Then, the problem of scattering by such 
grating also includes the previously considered case of a periodically modulated conductivity. 



In this case the funetions Mp*-* and M^||p_„ can be written as 



Ml*) = 



ikzh 



(-l)"«m||«^ 

> m||n 

,(-1)™k™,|„V2 _ (_;L)p-ng(-l)'"+'K„||„/V2' 



M, 



m| |p— n 



(300) 



(301) 



Contrary to the case of the sine profile, the Rayleigh method for the sawtooth profile 
is neither convergent for some values of ft-, for a given w, nor for a given h for all 
values of w. We have checked that partial convergence requires h < _D/10. The 
convergence over an energy range has to be checked case by case. 



In Fig. 35 we plot the absorvance, reflectance, and transmittance for a grating 
with a sawtooth profile where the convergence could be achieved. As in the case 



of the sinusoidal profile discussed in Sec. 9.4 we see the presence of resonances 
associated with the coupling of the impinging ER with the graphene SPPs. For 
Q = ^ the resonance is located at hu) « 11.6 meV. As before, the value can be 
predicted from: 



ei + £2 



-Ephcq. 



'p 



^^0.452 x 0.2|« 11.8 meV. 



(302) 
(303) 



At finite angles of incidence, the single peak seen at 9 = splits into two peaks 
because the ATR scan line intercepts the n — —1,1 bands at different energies. 
The prediction for the positions of the two peaks follows from the calculation of the 
band structure of the SPPs for the sawtooth profile. It is interesting to note that 
the peak splitting does not occur in the case of the sine profile. 
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Fig. 35. Absorbance, reflectance and transmittance for the sawtooth proflle at two angles of 
incidence, 6 = and 6 = 45°. The dashed lines refer to the non-corrugated limit for 6 = 0. Here 
h/D = 1/15 and F = 0.6 meV; other parameters as in Fig. |32| 

10. Graphene on a metallic grating 

10.1. Relation between the conductivity and the dielectric function 
of a 3D metal 

So far we have discussed the resonant couphng of ER to SPPs in graphene lying 
on a dielectric grating. In this section we shall consider the case where the under- 
lying grating is made of a conductor, be it a metal, a doped semiconductor, or a 
mesostructure li^SI 

Let us first derive a simple relation between the conductivity and the dielectric 
function of a 3D metal, based on Drude's model 1^^. The equation of motion of an 
electron in a metal subject to an oscillating electric field is 



X dx 



m,— + TOe7^ = -eE{uj)-"^* , (304) 



where rrie denotes the electron mass and —e is the electron charge. Substituting 



X = xoe , into Eq. (304) we obtain 



xo^ ^f"\. . (305) 



The polarization is defined as P{uj) = —erieXo (units of C/m^), where is the 
electronic density (per unit volume) of the gas. The electric displacement field is 
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defined as 



D{uj) ^ eoe{uj)E{ijj) = eoE{uj) + P{uj) 



Therefore, we have 



eoEi.) - J^Ihm. . (306) 



e(w) = 1 ( ■ (307) 



On the other hand, Drude's model for the conductivity reads I^^Zl 



^T3D = , . , . (308) 

me(7 - i(jj) 



Comparing Eqs. (307) and (308) we obtain 



e(w) = l + i— . (309) 

£01^ 

If the dielectric screening by core electrons of the metal atoms is taken into account, 



the unity in Eq. ( 309 ) must be replaced by a background dielectric constant > 1 
123. For instance, for gold £00 — 1.53 while for most semiconductors it is of the 
order of 10. 

Thus, when dealing with a 3D metal the conductivity enters the problem through 
the dielectric function. Often it can be assumed that 7 ^ w, leading to a purely real 
dielectric function, that is, the metal acts as a dispersive dielectric with a dielectric 
constant that can negative. In this limit, the dielectric function has the simple form 

e{uj) = e^-'^, (310) 

where uS^ = n^.e^ / (me^o) denotes the plasma frequency of the bulk conductor. Al- 
though it is, of course, measured in s^^, for convenience we shall use values in meV, 
which correspond to plasma energy hjjp. 



10.2. Surface plasmon polaritons at a dielectric-metal interface 

Before discussing the form of the spectrum of SPPs when graphene is placed on 
a bulk conductor, we analyze the simpler case of a dielectric-metal interface. We 
assume a system where a dielectric of constant permittivity £2 is in contact with a 
bulk metal of dielectric function e{uj). The dispersion relation can be obtained from 



Eq. (47) by putting a = and substituting ei by e(a;), namely 



Kl K2 

n,2 



0, (311) 



where k{ — q — e(aj)w /c and K2 = <? — £20; /c . As a result, for real values of 
Kl and K2, e(w) has to be real and negative, implying uj < ujp. Sometimes in the 
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literature the eigenvalue equation (311 1 appears in the form 

which can be easily obtained from Eq. (iBll). Note that, since the argument in 




the square root above must be positive, it is necessary that e{uj) < —£2 and uj < 
^p/V^oo + £2, which gives an upper bound for the SPP frequencies. 



We can solve Eq. (311 ) for u 



U! — 1 

2eoo 2 £2^00 



- -^Jlcu^e^ + c2g2(g^ + - ^e^e^c^q^OJ^ . (313) 
^^2^00 * 

The minus sign in front of the square root is necessary to guarantee that uj < 
+ £2- This relation takes simple forms in two limiting cases. When (7 — > 0, 

we have 



cV _ cV(£2_+eoo)' 
£2 4e^eoow2 



^2^^^ ^.^ VC2^. ; ■ (314) 

'Ic^t (,i2 ^ ' 



In this case, the dispersion relation of SPPs is always below the spectrum of the 
ER in the dielectric, a situation that does not occur for SPPs on graphene (recall 
Sec.|4]). 

For g — > cx) we obtain 

0.^ = ^^- , , ,22 ^ (315) 

+ £2 4£oo(,£oo ^" £2jC ^7 

that is, the dispersion relation is always below lOpj \/toa + £2- Clearly, the dispersion 
curve of the surface plasmon-polaritons at the surface of a 3D metal is quite different 
from that for a 2D conductor (see Sec. 141). 



10.3. Graphene on a metallic grating 

As known, for frequencies below the plasma frequency a bulk metal reflects all light 
impinging on it, a consequence of the negative value of its dielectric permittivityj^If 
the surface of the metal has a periodic corrugation, the energy of the reflected wave 
is distributed among the specular and the diffracted orders, as seen in Fig. [36] In this 
figure we depict three reflectance efficiencies, TZq, TZi, and TZ-i, for normal incidence 
of the impinging radiation. As expected for normal incidence, we have TZi = TZ_i. 
Since we are considering the limit of a dispersive metal with no absorption, there is 
no dissipation due to SPPs. The rich structure seen in TZq, TZi, and TZ-i spectra is 
due to the distribution of the energy through the different diffraction orders. The 



"Strictly speaking, this is true for e^c = 1, zero damping and neglecting interband or impurity- 
related optical transitions that may take place in the same spectral region. 
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Efficiencies ^ _^ for a metal sin-grating with no graphene 

D=1 nm; e,=1 ; to =0.5 eV -- normal incidence 







I 1 


1 ' 1 












— flat metal 
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flat metal 
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R, 
R, 
T„ 



Fig. 36. Reflectance and transmittance of a sinusoidal metaUic grating (without graphene). We 
depict three reflectance efliiciencies, TZq, TZi, and TZ—i. The specular transmittance To is also 
depicted as a dashed line; To 7^ only for u) > uip. The parameters arc; €1 — £00 = 1, Sojp = 0.5 eV, 
D = 10 Aim, h = D/IQ (top), and h = D/20 (bottom). 



interpretation of the figure is rather difficult because the problem is non-linear in the 
frequency lo. For energies below 0.12 eV (smaller than the plasma frequency taken 
equal to 0.5 meV) only the specular (T^-o) order is propagating. Since the absorp- 
tion of the metal was neglected the reflectance is equal to unity; this corresponds 
to frequencies uj/c < G. Above 0.12 eV the orders n ~ ±1 {Tl±i) also become 
propagating and we have Bragg diffraction (orders with |n| > 1 are evanescent; this 
corresponds to frequencies G < oj/c < 2G). In this case SPPs with the wavevectors 
±G, excited by the incoming light thanks to the presence of the grating, become 
radiative (non-evanescent) and the diffracted energy is distributed among the three 
orders n = 0,±1 in a non-trivial way. In this case we have TZq -\- TZi + TZ-i = 1, 
since there is no dissipation in the metal. Above the plasma frequency the grating 
becomes partially transparent. As the frequency of the incoming ER increases the 
dielectric function of the metal tends to Coo and To ^ 1, as seen in Fig. [36| 

If we deposit graphene on top of a corrugated metal surface, the system has two 
different regimes. For frequencies smaller than tUp/^/eoo + ^2 the system behaves 
essentially as the surface of a bulk conductive system. When u > cop/y/eao + ^2 the 
system behaves as a graphene sheet on a dispersive dielectric. The problem of the 
plasmon spectrum of a graphene sheet in the vicinity of a thick plasma-containing 
substrate was considered by Horing^iSSi^ -^^ho derived the dispersion relation of the 
surface plasmons of the system. 

In what follows, we assume that graphene is deposited on a metallic grating, as 
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illustrated in Fig. [37 




graphene 



Fig. 37. Graphene on a metal. The upper and lower media have dielectric permittivities £2 and 
ei = e(i^), respectively. 



For the metal we take a dielectric function corresponding to the limit of a disper- 

1. 



sive conductor, that is, given by Eq. (310). For simplicity, we shall assume e, 



We want to study the form of the dispersion relation of the surface plasmon- 



polaritons. In this case, the eigenvalue equation has the same form as Eq. (471 

K2e(w) 



1 



^162 



K2 
'^6062 



= 0. 



(316) 



Let us approximate the conductivity of graphene by its imaginary part only (the 
dispersive conductor limit). It allows for writing the eigenvalue equation as 



K2e(w) a AEp hcK2 
Kie2 €2 huj hiv 



0. 



(317) 



M=0.45 eV; r=0: t =2; e,=1 ; (o =0.04 eV 




q (1/nm) 



Fig. 38. SPP dispersion curves for graphene on a conductor with cOp = 40 meV and eoo = 1. The 
dashed black line is the light cone, q/ -y/c^ , and the horizontal line marks the value lu = ujp/ yj\ + £2- 
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Equation (317) gives the SPP dispersion relation for graphene placed on a metal. 

we 



When the conductivity of graphene vanishes, we recover Eq. (311). In Fig. 38 



plot the solution of Eq. (3171 for Up = 0.04 eV. The general trend of the solution 
is the following: when a; < Lo^j \J\ + 62, the SPP dispersion relation is that of the 
conductor underneath graphene; when lo > ujp/ \J\ + ti there is a change of the 
regime and the dispersion curve follows that of SPPs in graphene cladded by two 
dielectrics. It should be noted that, for a bulk conductor alone, the SPP dispersion 
relation lies below the value ujpj \J\ + £2 for all q (see previous section). We also note 
that the spectrum for the combined system { graphene+conductor } lies above that 
for the graphene alone (with ei = 1) because for w > we have < t(uj) < 1. 
Indeed, when uj > LUp and g 3> w/c, we can find an analytical expression for uj{q). 



In this regime, Eq. (317) reads: 



+?:^^ = o. (318) 



£2 eowe2 



Solving it in order of w gives 




^=\^+4f^^. (319) 



If Wp — > 0, we recover the dispersion relation of SPPs on a graphene sheet cladded 
between two media of relative permittivities ei = 1 and ti. 



In Fig. 39 we plot the reflectance, transmittance, and absorbance of graphene 
on a conductive sinusoidal grating. For the plasma frequency we assume a value of 
20 meV, corresponding to a doped semiconductor (for example n— GaAs, with 10^^ 
electrons per cm^); some types of mesostructures can also have very low plasma 
frequencies '^221 ^pj-^g choice of such a low plasma frequency simplifies the analysis, 
since, in this case, the only propagating order is the specular one. If there is no 



graphene on the grating (central panel of Fig. 39 ), we have total reflection for 
a; < Wp = Wp/ye^, as it should be. When the graphene is present, the reflectance is 
smaller than 1 for a; < Wp because of the absorption in graphene. Note, however, that 
if graphene were supported by a dielectric instead of the conductor, the reflectance 



would be much lower (top panel of Fig. 39 ). Above Wp, for this choice of parameters 



the reflectance and the transmittance show coupling of ER to the SPPs in graphene. 
The presence of the conductor beneath the graphene sheet shifts the position of the 
peak of the resonance toward higher energies because the its dielectric permittivity, 
although positive, is smaller than 1 (note that eoo = 1 in this example). 

As it can be seen in Fig. [39j the presence of graphene introduces a dip in the 
transmittance, at w > Wp. At the same time, an enhancement of the absorption of 
graphene is observed. It is also clear (bottom panel) that the reflectance drops to 
zero at a given frequency above 0.03 eV, which does not occur when graphene 
stands alone on a dielectric grating (top panel). Note that the metal does not 
absorb electromagnetic radiation, since in this model its dielectric constant is purely 



12, 2013 2:29 WSPC/INSTRUCTION FILE UMPByQ 



71 



o)p=0.02eV, E =ej=l,r=0.6 meV, ep=0.45 eV, D=10nni, h=D/10 




Fig. 39. Reflectance, transmittance, and absorbance of graphene on a conductive sinusoidal grat- 
ing. Top panel: free standing graphene with a sinusoidal profile. Central panel: sinusoidal metallic 
grating without graphene. Bottom panel: graphene on a sinusoidal metallic grating. The parame- 
ters are: eg = too = 1, T = 0.6 meV, ojp = 0.02 eV, Ep = 0.45 eV, D = 10 /^m, and h = D/10. 

real. The absorption occurs entirely within the graphene sheet. If the conductive 
substrate also absorbes ER, the analysis becomes more complex 



11. Summary 

At present, graphene plasmonics is an active field of research, both theoretical and 
experimental. The community has been addressing a broad range of topics in this 
field, some of which are: (i) coupling of ER to SPPs in graphene; (ii) use of nano- 
emitters for excitation of SPPs; (iii) enhancement of light absorption in graphene 
owing to SPPs; (iv) filters and polarizers exploring the propagation of SPPs along a 
graphene-covered surface; (v) launching and detection of SPPs in graphene, aiming 
at nanoplasmonic-based circuitry; (vi) exploring graphene double-layers (or mul- 
tilayers) where two (or more) different SPP branches emerge; (vii) ER switches 
and polarizers based on the ATR configuration; and (viii) formation of polaritonic 
crystals with band gaps controlled by gate voltage. 

In this work we focused on the SPP spectra of single and double layer graphene 
systems. We showed that the ATR configuration can work as an optical switch and 
that the double-layer system allows for moving the resonant energy towards higher 
energies. We showed further that, exploring the different response of graphene to 
TM and TE waves, one can control the state of polarization of the outgoing wave by 
changing the electron Fermi energy in graphene. It was demonstrated that SPP ex- 
citation can be achieved directly by illuminating a thin metallic stripe deposited on 
top of graphene, which can be viewed as a topological defect in the otherwise uniform 
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system, which enables to overcome the restriction imposed by the in-plane momen- 
tum conservation. We discussed the physics of graphene-based gratings where the 
underlying material is cither a dielectric or a conductor. We also showed that, even 
in the absence of a grating, an efficient ER-SPP coupling exists in graphene systems 
with periodically modulated conductivity, which can be induced by different mech- 
anisms (e.g., strain, doping, and gating). In particular, it was shown for an array 
of microribbons and for a continuous graphene sheet with a cosine-modulated con- 
ductivity. Both systems have the properties of a polaritonic crystal, although there 
is some difference in details between them. Another system from which one should 
expect qualitatively similar properties is graphene deposited on a periodically corru- 
gated surface (a grating) . The presented computational description of structures of 
this type in terms of Rayleigh-Fano expansions is demanding because the problem 
is poorly convergent. Even though, the cases of a sinusoidal and a sawtooth grating 
were considered and we showed that the ER-SPP coupling manifests itself by res- 
onances in the reflectance, transmittance, and absorbance spectra. When graphene 
is placed on a metallic surface with smooth periodic profile, the SPP dispersion is 
hybrid, showing the properties of a bulk metal at low wavemimbers and those of 
an isolated graphene sheet as large wavenumbers. Not surprisingly, the change of 
the regime takes place close to the plasma frequency of the metal underlying the 
graphene sheet. 
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